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ABSTRACT 


In  this  thesis  we  have  obtained  the  index  of  refraction  and  the  attenua¬ 
tion  coefficient  per  unit  time,  as  well  as  per  unit  length,  for  an  arbitrary 
medium  using  Maxwell's  wave  theory,  and  expressed  the  results  in  terms  of  the 
microscopic  currents  due  to  the  motion  of  the  particles  in  the  medium.  In 
Maxwell's  theory,  the  effect  of  the  medium  is  characterized  by  a  conductivity 
tensor  which  relates  the  macroscopic  current  to  the  electric  field.  The  con¬ 
ductivity  tensor  is  obtained  quite  generally  by  using  Kubo's  linear  response 
theory  in  terms  of  the  microscopic  currents. 

The  index  of  refraction  and  the  attenuation  coefficient  per  unit  time 
associated  with  the  decay  of  EM  waves  in  time  for  weakly  absorbing  media  has 
been  found  to  be  the  same  as  those  obtained  by  using  the  photon  transport 
theory.  However,  the  index  of  refraction  and  the  attenuation  coefficient  per 
unit  length  associated  with  the  decay  of  EM  waves  in  space  are  different 
from  those  obtained  by  the  transport  theory  even  for  weakly  absorbing  media. 

We  have  also  investigated  in  this  thesis  the  contribution  of  the  neutral 
atoms  to  the  absorption  of  photons  in  plasmas,  by  extending  the  Akcasu  and 
Wald's  work  on  the  absorption  due  to  the  inverse  bremsstrahlung  of  electrons 
in  the  field  of  neutral  atoms  to  higher  electron  temperatures  and  higher 
photon  energies,  and  formulated  this  problem  by  using  plane  waves  for  electron 
wave  function.  In  this  way  we  have  obtained  an  expression  for  the  absorption 
coefficient  per  unit  length  due  to  the  above  absorption  mechanism  in  terms 
of  the  elastic  and  inelastic  electron-atom  scattering  cross  sections  allowing 
the  atoms  initially  to  be  in  any  excited  state.  We  have  calculated  the  absorp¬ 
tion  coefficient  explicitly  for  hydrogen  atoms,  and  presented  the  results 
graphically  as  a  function  of  electron  temperature  and  radiation  frequency. 

Using  these  curves  and  the  conventional  formula  for  the  absorption  due  to  the 
photoionization  and  its  inverse,  we  have  computed  the  net  total  absorption 
due  to  the  neutral  atoms  numerically,  and  compared  our  results  to  the  absorp¬ 
tion  measured  by  Litvak  and  Edwards.  In  estimating  the  distribution  of  the 
neutral  atoms,  as  well  as  the  size  of  the  plasma  produced  by  the  laser  pulse 
in  their  experiment,  we  have  used  the  point  explosion  theory  with  spherical 
shock  wave.  The  agreement  between  the  calculated  and  measured  absorptions 
has  been  found  to  be  better  than  a  factor  of  10  and  in  fact  better  than  a 
factor  of  6  in  all,  but  one,  initial  gas  pressures  (The  observed  discrepancies 
may  be  attributed  mainly  to  the  use  of  the  radius  of  the  peak  luminous  volume, 
which  is  assumed  to  be  shperical,  as  the  actual  shock  wave  radius). 

In  the  absence  of  any  accurate  information  for  the  plasma  size,  and  of  an 
explosion  theory  which  takes  into  account  the  finite  initial  volume  of  the 
explosion  caused  by  the  laser  beam,  the  agreement  obtained  is  considered  as  a 
strong  evidence  for  the  importance  of  neutral  atoms  in  certain  absorption 


vii 


experiments  in  plasmas  over  ions,  because  the  absorption  calculated  by  con¬ 
sidering  the  electron  inverse  bremsstrahlung  in  the  field  of  ions  only  is 
about  100  times  less  than  the  observed  values. 


viii 


INTRODUCTION 


This  thesis  contains  primarily  an  investigation  of  low-energy  photon 
absorption  in  an  arbitrary  medium.  At  low  photon  energies,  the  pair  pro¬ 
duction  and  annihilation  processes  are  negligible  and  scattering  is  char¬ 
acterized  by  the  Thomson  cross  section  rather  than  Compton  cross  section. 

In  stellar  systems,  the  Thompson  scattering  may  be  important,  but  in  lab¬ 
oratory  plasmas  it  is  usually  negligible  because  of  its  small  cross  section. 

Thus  our  analysis  in  this  thesis  is  restricted  to  the  systems  in  which  the  scat¬ 
tering  is  negligible,  and  the  net  photon  absorption  is  primarily  due  to  brems- 
strahlung  and  inverse  bremsstrahlung,  ionization  and  recombination,  and 
excitation  and  deexcitation. 

The  equations  of  radiative  transfer  for  dispersive  and  nondispersive 


media  are  usually  developed  by  phenomenological  consideration. 


(1-10 


systematic,  self-contained  derivation  of  a  photon  transport  equation  for 

nondispersive  media  from  fundamental  consideration  was  made  in  1961  by 

(5) 

Osborn  and  Kievans.  In  their  derivation,  a  photon  distribution  func¬ 
tion  in  analogy  with  the  quantum  mechanical  distribution  function  for 
(6) 

particles  was  introduced.  The  equation  they  obtained  can  be  reduced  to 
the  conventional  radiative  transfer  equation  obtained  phenomenologically 
for  nondispersive  media. 

(7) 

A  year  later,  they  extended  their  theory  to  dispersive  media  by 

(8) 

making  use  of  the  concept  of  "dressed  photon"  first  introduced  by  Mead. 


1 


2 


A  "dressed  photon"  has  a  different  frequency  in  the  medium  than  the 
free  space  value  ck,  but  it  has  the  same  wave  length  in  the  medium  and 
in  free  space.  With  the  "dressed  photon"  technique  they  derived,  in  the 
framework  of  the  first  order  perturbation  theory  a  photon  transport  equa¬ 
tion  for  dispersive  media.  This  equation  differs  from  the  radiative  trans- 

(3  4) 

fer  equation  obtained  phenomenologically  in  dispersive  media  by  others.  * 

(9) 

In  order  to  compare  these  two  theories,  Wald  performed  an  experiment  in 
1966  in  which  he  measured  the  absorption  of  microwave  radiation  in  slightly 
ionized  helium.  Although  a  better  agreement  is  obtained  by  the  photon 
transport  equation  than  the  radiative  transfer  equation  in  his  measurement, 
conclusive  evidence  of  the  validity  of  the  photon  transport  equation  or 
the  invalidity  of  the  radiative  transfer  equation  cannot  be  inferred  from 
his  measurement  because  the  refractive  defocussing  effects  are  not  neglibible. 

Since  the  net  photon  absorption  and  the  refractive  index  in  dispersive 
media  can  be  easily  deduced  from  the  photon  transport  equation,  the  validity 
of  this  equation  can  be  tested  if  the  net  photon  absorption  and  the  refrac¬ 
tive  index  can  be  obtained  independently  by  an  entirely  different  approach. 
Such  an  approach  to  the  calculation  of  the  net  photon  absorption  and  the 
refractive  index  can  be  achieved  by  using  Kubo's  theory^0 ^  for  electric 
conductivity  and  the  Maxwell  equations.  In  this  approach,  the  absorption 
coefficient  and  the  index  of  refraction  are  expressed  in  terms  of  the  mi¬ 
croscopic  current  due  to  motion  of  all  particles  in  the  medium.  To  facil¬ 
itate  the  comparison  we  obtain  in  section  II- 1  the  net  photon  absorption 
and  the  refractive  index  using  the  photon  transport  equation  also  in  terms 
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of  the  microscopic  current.  The  expected  value  of  electric  current  is  ob¬ 
tained  from  Kubo's  theory  in  section  II-2.  Using  the  expected  value  of  the 
current  in  the  Maxwell  equation  (II-40),  a  dispersion  relation  between  the 
wave  vector  k,  the  frequency  cu,  and  the  electric  conductivity  can  be  ob¬ 
tained.  Two  different  sets  of  results  for  photon  absorption  and  refractive 
index  in  weakly  absorbing  media  are  obtained  by  considering  the  damping  of 
the  electromagnetic  wave  in  time  and  in  space.  A  comparison  of  the  results 
obtained  from  these  two  theories  is  presented  in  section  II-3. 

In  order  to  display  the  various  mechanisms  contributing  to  the  photon 
absorption  and  also  to  estimate  the  order  of  magnitude  of  the  various  con¬ 
tributions,  wtf  use  a  convenient  and  simple  representation  for  the  particle 
system  in  Chapter  III.  Second  quantization  is  used  to  express  the  various 
potentials  between  particles,  as  well  as  the  interactions  between  particle 
and  radiation,  in  terms  of  the  particle  and  radiation  creation  and  destruc¬ 
tion  operators.  Starting  from  the  golden  rule,  one  can  obtain,  with  several 
approximations  discussed  in  section  II-2,  a  simple  expression  which  dis¬ 
plays  the  various  mechanisms  contributing  to  the  photon  absorption,  such 
as  the  bound-bound  transition,  bremsstrahlung  of  electrons  in  the  fields  of 
the  neutrals  and  ions,  the  induced  dipole  transitions,  etc.  At  the  end  of  this 
chapter  we  give  a  simple  and  crude  investigation  of  the  variation  of  photon  ab¬ 
sorption  with  time  after  the  formation  of  a  plasma.  This  investigation  is  mo¬ 
tivated  by  the  absorption  measurements,  being  performed  now  in  The  University 
of  Michigan,  in  which  a  continuous  He-Ne  laser  beam  (6328A)  is  incident  on  the 
decaying  plasma  produced  by  exploding  lithium  wire.  The  validity  of  this 


k 


simple  investigation  could  not  be  verified  because  the  measurements  have  not 
reached  the  final  stage  yet. 

In  Chapter  IV,  we  will  be  concerned  with  the  radiation  absorption  only 
due  to  inverse  bremsstrahlung  of  electrons  moving  in  the  field  of  neutral 
atoms  at  high,  gas  temperatures  (~  20  eV).  The  same  problem  at  low  gas  tem¬ 
peratures  (~  1  eV  or  less)  was  investigated  in  i960  by  Firsov  and  Chibisov^1^ 

(12) 

and  recently  extended  by  Akcasu  and  Wald.  At  low  temperatures,  the 

electron  energies  are  insufficient  to  excite  an  atom  from  its  ground  state 
to  an  excited  state.  Assuming  that  all  the  atoms  in  the  system  are  initially 
and  finally  in  the  ground  state,  they  calculated  the  various  absorption  con¬ 
tributions  due  to  neutral-inverse  bremsstrahlung,  induced  dipole  transition, 
and  exchange  and  interference  effects  by  partial  wave  method  and  found  that 
the  last  three  contributions  for  low  temperature  system  are  negligible  as 
compared  to  the  first  one. 

An  experiment  measuring  the  absorption  coefficient  for  the  ruby 
laser  beam  (69U3X)  in  a  hydrogen  plasma  produced  by  a  giant  pulsed  laser 
beam  was  carried  out  in  1966  by  Litvak  and  Edwards.  Their  calculated 

absorption  coefficient  obtained  by  considering  photoionization  and  inverse 
bremsstrahlung  of  electrons  in  the  field  of  ions  is  two  orders  of  magnitude 
less  than  their  measured  result.  Chapter  IV  is  motivated  primarily  by  this 
large  discrepancy. 

The  temperature  of  the  plasma  in  Litvak  and  Edwards  experiment  is  high 
(9  ~  20  eV).  Most  of  the  hydrogen  atoms  are  found  in  excited  states  and 
the  energies  of  electrons  are  sufficient  to  excite  the  atoms  from  a  level 
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to  a  higher  level.  The  electron-atom  scattering  cross  section  increases  when 
the  atom  is  in  an  excited  state.  As  seen  from  Akcasu  and  Wald's  work,  this 
cross  section  enters  in  the  expression  for  radiation  absorption.  Since  at  high 
temperatures,  an  appreciable  number  of  hydrogen  atoms  are  in  upper  levels  and 
the  above  cross  section  at  these  levels  is  high,  one  may  expect  the  neutral 
bremsstrahlung  to  be  a  dominant  process  contributing  to  the  radiation  absorption 
in  Litvak  and  Edwards'  experiment.  This  is  one  reason  for  extending  in  Chapter 
IV  the  calculations  by  Akcasu  and  Wald  to  high  temperatures.  The  other  reason 
is  that  the  total  number  density  of  hydrogen  atoms  in  different  states,  deter¬ 
mined  in  the  case  of  Litvak  and  Edwards'  experiment  by  the  initial  gas  con¬ 
ditions  with  the  assumption  of  ideal  gas,  is  ten  or  more  times  the  electron  (or 
ion)  density  depending  upon  the  initial  gas  pressure. 

In  the  formulation  of  the  problem,  we  use  the  second  order  perturba¬ 
tion  theory  in  which  electron  states  are  represented  by  plane  waves.  In 
this  approximation,  the  electron-hydrogen  cross  section  is  calculated  in  the 
first  Born  approximation.  Although  a  more  accurate  result  is  expected  by 
using  partial  waves ,  the  use  of  free  electron  wave  function  makes  the  prob¬ 
lem  more  manageable. 

If  all  the  atoms  in  the  system  are  initially  and  finally  in  the  ground 
state,  as  assumed  by  Akcasu  and  Wald,  with  the  assumption  that  the  cross 
section  involved  is  slowly  varying  up  to  the  incident  energy  of  electrons, 
we  obtain  the  same  result  as  obtained  by  Akcasu  and  Wald  through  the  partial 
wave  method  (section  IV-2). 


An  application  of  the  above  theoretical  results  to  hydrogen  plasma  is 
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presented  in  sections  IV-3,  IV-4,  and  IV-5-  In  order  to  explain  the  measured 
absorptions  in  Litvak  and  Edwards  experiment  by  considering  the  photoionization 
process  and  the  inverse  bremsstrahlung  of  electrons  in  the  field  of  neutral 
atoms,  we  use  in  section  IV-7  the  point  explosion  theory  in  estimating  the 
number  density  of  neutral  atoms  in  the  plasma  produced  by  the  laser  beam.  The 
comparison  between  our  calculated  absorptions  and  the  measured  results  are  dis¬ 
cussed  in  section  IV-8. 


CHAPTER  II 


PHOTON  ABSORPTION  AND  REFRACTIVE  INDEX 

In  this  chapter  we  shall  present  two  different  approaches  to  the  cal¬ 
culation  of  photon  absorption  and  refractive  index  in  an  arbitrary  particle 
system.  The  first  approach  is  based  on  photon  transport  theory  developed 
by  Osborn  and  Kievans ^>7)  secon(j  one^  on  Kubo's  theory of  electric 
conductivity.  The  comparison  between  the  results  from  these  two  theories 
will  be  made  at  the  end  of  this  chapter. 

1.  PHOTON  TRANSPORT  THEORY 

(5 ) 

In  1961,  Osborn  and  Kievans  used  first  order  perturbation  theory  in 

(7) 

developing  photon  transport  for  nondispersive  media,  and  a  year  later  they 

extended  the  theory  to  dispersive  media  by  making  use  of  the  concept  of 

(8 ) 

"dressed  photon"  first  introduced  by  Mead.  In  this  section  we  shall  use 
their  results  for  dispersive  media  to  express  the  photon  absorption  and  the 
refractive  index  in  terms  of  the  microscopic  current  due  to  motion  of  all 
particles  in  the  medium.  This  requires,  in  the  first  place,  a  description 
of  the  hamiltonian  of  the  particle  system,  the  radiation  field,  and  the 
interaction  between  these  two,  as  well  as  the  introduction  of  the  concept 
of  "dressed  photon." 

Consider  a  box  of  volume  L3  in  the  particle  system  under  consideration. 
The  hamiltonian  for  the  particles  in  this  volume  element  interacting  with  a 
radiation  field  can  be  written  in  the  nonrelativistic  theory  as 
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with 


d  ^  ^  l-oj  -  “ 

h  =  hR  +  ?  5 - I 


(£oj  - 

2moJ 


+  V0  +  V'  =  H0  +  H 


R  P 
Hq  *  H  +  H 

hPR  =  hPR1  +  hPR2 


(II-1) 


(II-2) 

(II-3) 


where 


+  VJ  +  V' 


H  =  Z  I 


H  =  "  c5  2m0jC  ) 

E®-AA2(rfl1) 

oj  2majC-  -  ~0J 


(II-4) 


(II-5) 


(II-6) 


In  the  above  equations,  and  A(r)  are,  respectively,  the  hamiltonian 
and  the  vector  potential  of  the  radiation  field.  The  symbols  m^,  e0j. 
pQj ,  and  rQj  denote  the  mass,  charge,  momentum,  and  position  of  the  j-th 
particle  in  the  e-th  molecule.  Here,  we  use  the  term  "molecule"  in  a  gen¬ 
eral  sense  to  refer  to  any  aggregate  of  particles  bound  together.  The  num¬ 
ber  of  constituent  particles  in  a  molecule  is  arbitrary.  It  proves  to  be 
convenient  to  regard  even  an  electron  as  a  simple  molecule  as  defined  above. 
VQ  is  the  potential  between  the  particles  in  the  e-th  molecule  and  V'  is  the 
potential  between  the  molecules  in  l?. 

Interactions  between  particles  in  different  boxes  via  long-range  coulomb 
forces  constitute  a  small  effect  upon  photon-particle  interactions  within  a 
box  and  we  may  expect  that  neglect  of  this  effect  produces  negligible  error. 
Then  the  photo-particle  interaction  of  the  system  can  be  approximated  as  the 


sum  of  the  interactions  in  each  box. 


mmrvptmm 


The  customary  procedure^*1  ^  requires  that  the  (transverse)  radiation 
field  be  periodic  at  the  boundaries  of  a  normalization  box  whose  volume  is 
assumed  here  to  be  1? .  The  vector  potential  of  the  field  at  a  point  r 
can  then  be  represented  by  operators  in  the  Schr8dinger  picture  as 


i(i) 


Br  e'1- 


(II-7) 


and  the  hamiltonian  of  the  radiation  field  turns  to  be 


H8  =  ^  ^  ck  {a^(k)a^(k)+Q^(k)a^(k)} 


(II-8) 


T  ° 

where  a^(k)  and  0^(k)  are  creation  and  destruction  operators  for  photons 
of  momentum  ftk  and  polarization  \  in  free  space.  £  (k)  is  a  unit  polarization 

“  A. 

vector.  The  creation  and  destruction  operators  satisfy  the  commutation  re¬ 
lations 

[c^(k),  C^,(k' )]  =  [o^(k),  0^,(k’)]  =  o 


[o^(k),  C^,(k'  )]  =  6^,6. 


kk' 


(II-9) 


HMSSOa 


Since  the  radiation  is  in  constant  interaction  with  the  medium,  Mead 
introduced  the  concept  of  "dressed  photon"  by  associating  photons  in  the 
medium  with  a  different  frequency  than  the  free  space  value  ck,  keeping 
the  wavelengths  in  medium  and  in  free  space  the  same.  In  doing  this,  he 


used  a  different  expansion  of  A(r)  due  originally  to  Bohm  and  Pines 


(15) 


and  showed  that  the  creation  and  destruction  operators  a  (k)  and  a  (k)  for 

K  K 

creating  and  destructing  photons  of  frequency  are  different  from  the  free 
space  operators  Ci^(k)  and  Q^(k).  They  are  related  through  the  relationship 
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\<*>  ■  i© 2  sr>\t(-k» 


(II-10) 


C^(k)  = - f  ((ck-KJok)a^(k)+(ck-a)k)a^(-k)} 

2(c%ck)^  (II- 11) 

Q^(k)  = - ((ck+aik)a^t(k)+(ck-oDk)a^(-k))  . 

2(%ck)^ 

It  is  clear  that  ( 1c )  and  a^(k)  become  identical  when  =  ck.  It  can  be 
verified  easily  that  a^(k)  ant^  ^  satisfy  the  same  commutation  relations 
as  0^(k)  and  a^(k)  even  with  ^  ck. 

The  substitution  of  eqs.  (11-11)  into  eqs.  (II-6),  (II-7),  and  (II-8), 

PR2  R  t 

i.e,,  expressing  H  ,  A(r),  and  n  in  terms  of  the  operators  a  (k)  and 

A.  “ 

ax(k),  gives 


"lL“k 

+  H^1  + 


(II- 12) 


(II- 15) 


HpE2  .  HoPR2  +  Hra2  +  HsPR2 


(ii-iM 


where 


to 

^  *4 T  (\(i)\(Js)*\(ls)\(£)) 

c2l&,  2 

n”  *4»“^T’  ,ak(->\(k)t\(iE)ak(-)) 

c^k ^ 

hR2  *4?  -zr 


(II-15) 


(11-16) 


(H-17) 


ct3>»  k  5 

"  L  moi<\ 


aj  t 


(a*(k)a^(k)+a^(k)a^(k)}  (11-18) 
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H 


PR2 


irfis 


sL 


aj\k  .3 


L  m 


A 


(a^(k)a^(-k)+a^(-k)a^(k) )  ( 11-19) 


«2PR2  -  Jr  k 


Jtft  ed  -i(k-k' )  'T 

__2i  ® -  ■  'i  a"  (a  (kR  (k)+a  (-k)e(-k)} 

l\j  (a^.Ja  *•  ^  ^ 


’  ( Sv  ( -k')iv  ( -k1) +6^  ( k1)  ( k1) } 


(11-20) 


i 

Jf'k' 


means 


\k¥\.'k' 


Then  the  hamiltonian  for  the  whole  system  (particles  in  I ?  plus  radiation) 
can  be  written,  in  terms  of  the  operators  a^(k)  and  a^(k),  as  H=H!+H*  with 
H0  =  HR0+HP  and  HI=HPR1+(HR1+H0PR2)+(HR2+H1PR2)+H2PR2  where 


Une 


flL 


b  \k)o3  T3„  ...  I  v  \ 


c2k^1  t  t 

*  "" sr  K  <JsK  (JsK  OsK  (is 


**  *  "i1*2  *  !&<  J 


L'V-“k  J 

We?.  c2kS»f 


x(k)+a^(k)a^(k))  (11-21) 


2JL_  + 


h  \Tt\a3  .3 

L \a 


(avt(k)at(-k)+a  (k)a (-k)}  (11-22) 

I  K  K  K  \ 


,(7) 


In  obtaining  the  photon  transport  equation  and  refractive  index,  Kievans' 
formulated  the  problem  in  the  representation  in  which  the  particle  state  |  n>  and 
the  photon  state  I  q>  satisfy 


HP|n>  =  Ejn>  HR°|q>  =  E^|q> 


i.e., 


(II-23) 


H’jnq>  =  E  |nn>,  E  =  E  +E 
o  nr)  nq  n  q 

With  HP  regarded  as  the  perturbed  hamiltonian.  The  obtained  photon  transport 


equation  is  given  by 
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~  ~t~  "  *  £'Vvfx(r,k,t)  -nJn.n,  (^t— )-\(->)  (H-*) 

where  f  (_r,_k,  t)  is  the  expected  number  of  photons  with  momentum  "hk  and 
\ 

3 

polarization  \  in  the  volume  L  located  at  the  point  r.  The  photons  are 
moving  with  the  velocity  vfi  in  the  direction  of  the  unit  vector  fi.  v  is 
different  for  dispersive  media  from  the  speed  of  light  c.  |  n  '  t) ’ ,  n-rj  is  the 
transition  probability  from  the  initial  state  |nq>  to  the  final  state 


I n ’  -rj ’  >  which  is  given,  in  first  order  perturbation  theory,  by 

Tn,Tl>T1--f  |<n,r1,|HI|nn>|2&(En,Ti,.EnTi).  (H-25) 

r^(k)  is  the  occupation  number  of  photons  with -ft  k  and  and  D  is  the  den¬ 
sity  operator  of  the  whole  system  (particles  and  radiation). 

In  addition  to  the  photon  transport  equation,  the  refractive  index  of  the 
medium  can  also  be  obtained  by  letting 


F  D  (S  -S  )  =  o 

nq  nq,nq  nq  no' 


(11-26) 


where  SnT  is  the  shift  of  the  energy  level  E^  for  the  state  |nq>  and  given 


snn=<„,iHW  ♦  i  p  i^vihWi8 

1  nV/nr,  Enq-En'q' 


(H-27) 


where  P  indicates  the  principal  value.  Sno  is  a  shift  of  the  self-energy  of 
the  medium  when  no  photons  are  present  and  can  be  obtained  by  letting  q=o 
in  eq.  (11-27) • 

With  the  above  results,  we  shall  express  the  photon  absorption  coeffi¬ 
cient  and  the  refractive  index  of  the  medium  in  terms  of  the  microscopic 

3 

current  due  to  motion  of  all  particles  in  L  in  the  following  two  sections. 
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a.  Photon  Absorption  Coefficient 

The  processes  involved  in  obtaining  the  photon  absorption  coefficient 

are  single  photon  emission  and  absorption.  The  only  contribution  to  these 

processes  comes  from  H  because  n  ,  HQ  ,  Ji  ,  ,  and  Hg  are 

bilinear  in  photon  operators .  Define  the  current  operator  due  to  motion 

3 

of  all  the  particles  in  L  as 


(11-28) 


PR!  *  A(r) 

HPR  =  -  /d5r  J(r)  * 

where  the  integration  takes  over  the  volume  l? .  From  eqs.  (11-23)  and 


(II-29) 


the  definition  of  Dirac  &- function 


5(I)  -  £  C '  11  d-> 

the  substitution  of  eqs.  (11-29)  and  (11-25)  into  eq.  (11-24)  gives 


at  (i>is,t) 


+  n  •  Wf^(r,k,t)  =  Q^(r,k,t) 


where 


P  _  _  A(t) 

ft  (r,k,t)  *  Z  /d^r  /d  r*  Jdt*  {  Z  <  n  |  J(r)  |  rf>  *  <t)|=-=-|  tj+  1> 

*  nTl  *2  ^  n’  c 

A(r *t 1 )  A(r) 

BKq+ll””” - h>*  <n' |  J(£»  t 1  )|n> -2J(<n|  J(r)  |  n”  >  *  <tiI=T=”I  n- 1  > 


(II-30) 


A  |  j*  *  ^  *  j 

■  <T)-l|— - |tj>  •  <n"|j(r;t'  )|n>  ) 


(11-31) 


hpt  -i.Hpt 
j(r,t)  ■  J(r)e  * 


(II-32) 


11+ 


i^ot  .iHROt 

A(r,t)  =  A(r)ett  (II-53) 

being  the  Heisenberg  operators.  The  fir.t  sum  in  the  bracket  of  eq.  (II-3I) 
comes  from  photon  emission  by  letting  q'  =  q+lin  eq.  (11-24).  The  second 
term  comes  from  photon  absorption  by  letting  q'  =  q-l. 

Evaluating  the  matrix  elements  of  A(r)  and  A(_r,t)  in  Eq.  (H-31)  by 
use  of  eqs.  (11-12)  and  (11-33),  and  suming  the  intermediate  states  |n’> 
and  |n">,  ore  obtains 
2«D , 


0  (r,k,t)  =  Z  "-r}T];-nT]  /d5r  /d3r'  /dt '  <  n |  J  (r)J  (r,’t ' )  |  n  > 
^  -  nrl  ^  ^ 


b  ^e 


where 


io)  t'  ik-(r-r') 


^(r,k)+l  ^(r,'k)+l- 


-io)  t'  -ik  •  ( r-r  '  ) 


J^(r,t)  =  J(r,t)’e^(k). 

At  tnis  stage  we  introduce  several  approximations  which  enable  us  to 

reduce  Q^(r,k,t)  to  a  simple  form.  The  first  approximation  is  to  replace 
PR  P 

^nq,nq  by  D  nn  D  qq*  ^  the  density  operator  of  the  medium  only  and 
P  _pj{P  1 

given  by  D  =e  /Tre  with  P=^  being  the  reciprocal  of  the  medium  tern- 

perature  in  the  units  of  energy.  D  is  the  density  operator  of  the  radia- 

R  -  PH^®  -  BH^ 

tion  field  and  given  by  D  =e  /Tre  .  This  approximation  implies 
that  the  particle  system  and  the  radiation  field  are  initially  statistically 
independent  and  permits  us  to  perform  the  statistical  averages  over  the 
particle  and  photon  states  separately.  The  second  approximation  is  equiv¬ 
alent  to  replacing  the  average  of  a  function  by  the  function  of  the  average 


in  performing  the  statistical  average  of  the  factor  containing  the  photon 
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occupation  number,  i.e., 


Z  Dqq  Jl+i^(r,k)  '  Jl+T^(r;k),te  j  l+fx(r,k)'  J  l+f^r.'k), 

With  these  two  approximations,  t)  reduces  to 


(iio  t  *  ,  ,  ,  . 

e  k  eik  L~L  J  l+fx(r,k)'  J  l+f^(r.'k) 

r. 


t') 

We  recall  that  integrations  over  r  and  r'  are  extended  to  the  same  box  of 
volume  l? .  If  we  further  assume  that  f  (r,k)  is  slowly  varying  over  the  box 
volume  then  we  can  take  the  product  of  the  square  roots  outside  the  integral. 
Interchanging  r  and  r'  and  letting  t'=-t'  in  the  second  term,  one  finally 
obtains,  by  using  the  property  TrABC=TrBCA, 


-iV’  -lk-(r-r') 
-e  e  - - 


f*(l>k)fx(l>k) 


TrDP  Jx(r) 


J|r  +  £*Dvfx(r,k,t)  =  -o^(k,u>k)f^(r,k,t)  +  E^Ck,^)  (11-35) 


where 


<\(k,<\)  =  I  d3r  Jdr'ldt'e  k  TrDP[Jx(r,’t '  ),  J^r )  ]  (ll-3^a ) 

Lft\ 


2n 


L3tok 


/d5r  /d5r' 


7dt ' 


iix>  t 1  ,  . 

;  k  eik-(r-r  )  TrDp  j  (£)j  (r;t<)  (11-35) 


are  respectively  the  net  absorption  coefficient  and  the  spontaneous  emission 
coefficient  per  unit  time  for  photons  of  momentum  tik  and  polarization  In 

p 

an  infinite  homogeneous  medium,  TrD  [J  (r,'t' ),  J  (r))  can  depend  only  on  the 
difference  of  the  position-  r-r'.  In  a  large  finite  system,  this  translational 
invariance  will  be  approximately  true  in  regions  away  from  boundaries.  One  of 
the  integrations  over  the  position  can  be  performed.  Then 
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^“k)  7dfe  k  TrD?  [J^  ( r,’  t '  ),  ( r )  ]  (11-3*0 

In  view  of  eq .  ( 11-33)  it  is  obvious  that  the  photon  absorption  coefficient 


per  unit  length  is  given  by 


(11-36) 


In  a  later  section  we  shall  compare  eqs  ( 11-3*0  and  (II-36)  to  the 
expressions  for  photon  absorption  obtained  by  Kubo's  theory.  The  superscripts 
s  and  t  over  the  absorption  coefficients  per  unit  length  and  per  unit  time, 
respectively,  are  introduced  here  to  facilitate  this  comparison. 


b.  Index  of  Refraction 

We  mentioned  before  that  the  index  of  refraction  in  dispersive  media 
can  be  obtained  by  letting  DnT)>tlT,(snT)-Sno  )=°-  From  ed-  (lI-2 7), 


nj]  Dnn,np  (Snr|”Sn0  )  ^ 


Dnn,m[<"'llHl|"1>-<no|HI|no>  + J  P J<nEl'^ 

L  1  '  nq  n  r. 


<n ' n'  H  no  > 


p  1  11  '1  Lii  1 
n'q'/no  ^ho-En' 


(II-37) 


Recalling  H  and  Eqs  (II-5),  (11-20),  (II-21),  and  (11-22),  it  is  readily 


established  that 


<n  n  ( H1 1  n  r|  >  =  <„ti|  H^1+HC|PR2|  „q> 

-  1  J  urn  .CD. 

^  SJ  k 


IttN  e2.  c2\3a? 

s  s.i  ^ 


(2n(k)+l)  (11-38) 


.  1.J1  1  „R1  PR2,  .  r  ft  Tr  hNses,1  c2k’U\ 

<no  H  no  >  =  <no  H  +  H_  no  >  =  .  —  T  - “■  +  - 

0  %*^  si  3  CD. 

L  ms<A  J 

where  N„  is  the  number  of  the  molecules  of  kind  S  in  1?  and  e  .  and  m  . 

s  sj  sj 


(II-39) 
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are  the  charge  and  mass  of  the  jth  particle  in  a  molecule  of  kind  S. 

Note  that  <n'  q  |H^|nq>  =  <n'o  jH*|no>  =  o  for  n'/n,  and  <nr)'  |H^|nrj>  =  o 
for  tj'  =  r)±l  and  for  T)  =  o,  ri'  =  1,  then  one  obtains 

£  l<n 1  n*  I  H1 1  nn  >  |  2  _  y  |<n'n|HIlnn>  |2  +  y  l<nn*  |  H1 1  nn  >  |2 


n'r]'j/nrj  ^nri'En'^' 


n'/n 


^n“^n ' 


ErEn' 


+  y  l<n,n'|HI|nn>  l2  =  £  |<nV  |HI|nr)>  |2 

n'^n  Enn-En '  rj  *  nVn  Enrj-En'ri' 

TlV'H  ’lA 


(11-10) 


and 


£  |<n 1 n1 I H1 1 no  >  | 2  ^  [<n 1 n’ IH1 1 no  >  | 2 

n't)  '/no  Eno-En'^'  n'*n  En0-En'^' 

n'fo 


(II-ll) 


Using  the  property  of  Dirac  6-function  f(E)=  /  f (x )&(x-E)dx,  eqs .  (II-IO) 

-00 

and  (II-ll)  become 


z  l<n-.n.'l1|1ln'i>  1 2  g  z  /to. 

n'  rj'/nTj  EnTi_En’n'  n'fn-c 

*i?n 


.  l<n' n' jHFR1|nrt>  l\f  EnirEnV 


and 


li-j)' 


tPR1, 


6(u>' 


) 


r  lno>l  =  z  Jin-  |  l"°>]  6(a)._  ^  -'-I'l'  )  . 

n'qVno  ^no'^n'rj'  n'^n -»  *0)  -fi 

The  substitution  of  eq.  (II-29)  and  the  use  of  eq.  (II-3O)  gives,  after  a 
straightforward  manipulation  (see  Appendix  A), 


v  |<n  '  q'  |  H^l  nq  >  [  2  <.  |<n  '  q1  |H^[  no  >  [ ' 

n' q'^nq  ^nq"^n'q'  n' T)  /no  ^no"^n'q' 


\(E)  3  3  oo  oo  ik  •  (  r-r '  ) 


l\ 


-00  -00  CD 


A 


i<n|  Jv  (r'  )JX  (r,  -t '  )-J^(r  )J^(r,' t '  )|n>  . 


(11-1*2) 
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By  taking  the  statistical  average  over  the  photon  state  |  rj  >  a s  in  the  case 
of  absorption  coefficient,  using  TrABC=TrBCA  and  substituting  eqs.  (II-38), 

(II-39),  and  (11-42)  into  eq.  (lI-37)>  one  obtains 

InN  e^ 

r  s  s.1  2.  2  2 

h  k  7 —  c  * 

E-^fUr.lUO3  L  *..1 _ 


-'<£ ,TrDP  -»• 


^  2  - zj! — 

\k  2 

v.  k 


2  *  x  <»  00 

~ — —  /d  r/d  r ' /dt  P/do> ' 

J  2  — OC  —OO 

iTttU 

k 


Then  the  refractive  index  for  dispersive  media  is  obtained,  in  terms  of  the 

3 

microscopic  current  due  to  motion  f  all  the  particles  in  L  ,  as 


2 

no  =  (o>K)  =  1+ 


Unb^Qc,^) 


(11-^3) 


where 


iuPt '  2 

b  (k,o).  )  =  fdr'Jdt '  P  /do)'  ^ }TrDP[J  (r,'t’ ),  J  (r)]  (11-44) 

K'-*  k  2nfMy.  J  -s>  -*  cn'-u^  \  \  ~ 


2  Nsfsi 

q 


(II-U5) 


and  the  property  that  TrD^[J  (r,'t'),J  (r)J  can  depend  only  on  r-r'  for  an 
infinite  homogeneous  medium  has  been  used.  In  section  II-3,  we  shall  compare 
eqs.  (II-3I),  (II-3C),  and  ( 1 1- U3 )  with  the  results  obtained  from  Kubo's 


theory. 


2.  KUBO'S  THEORY 

(10) 

In  1957,  Kubo  developed  the  theory  of  linear  response  of  a  medium 
to  an  external  field  (i.e.,  electromagnetic  waves)  acting  on  the  medium. 

He  showed  that  the  response  can  be  described  by  electric  conductivity  tensor 
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of  the  medium.  In  1966  Dong^^  applied  the  theory  to  fully-ionized  plasmas. 

In  this  section  we  shall  use  Kubo's  theory  to  obtain  the  photon  (or  radiation) 
absorption  coefficient  and  the  refractive  index  for  an  arbitrary,  electrically 
neutral  medium  which  may  be  a  neutral  gas,  or  partially  or  full-ionized  plasma. 
From  Maxwell's  equations,  one  can  obtain 

2 

^E(r,t)  -~~E(r,t)  =  “  ^  J6(r,t)  (11-16) 

c  St  c 

by  assuming  that  the  macroscopic  charge  density  of  the  medium  is  zero,  i.e., 

V*E(r,t)  =  o. 

Thus  E(r, t)  is  transverse.  We  used  the  superscript  e  over  J(r, t)  in  eq. 

(II-i+6)  to  denote  that  it  is  the  expected  value  of  the  current  operator 
J(r, t).  By  using  the  gauge  in  which  the  scalar  potential  vanishes,  the 
electric  field  E(r, t)  can  be  described  by  a  vector  potential  A(r, t) 
through  the  relationship 


E(r,t)  =  -  -  —  A(r,t).  (11-1*7) 

Assume  that  the  external  field  turns  on  at  t=-«.  Before  the  field  is 
imposed  on  the  medium,  we  have 

[DP,  HP]  =  o. 

At  time  t,  the  interaction  hamiltonian  between  the  system  and  the  applied 
external  field  can  be  written  as 


=  H1(t)+H2(t) 


(11-1*8) 
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where 

*  A(  r,  t ) 

H  (t)=-/d5r  J(r ) •  = - 

i  c 


H2(t)=-  £  /d5r  JA(r,t)* 


A(r,t) 

c 


(11-49) 

(11-50) 

(H-51) 


With  J(j?)  defined  as  eq.  (11-28).  When  the  medium  is  acted  upon  by  the  field 
A(r,t),  the  total  current  operator  at  the  time  t  becomes 

J(r,t)  =  J(r)+JA(r,t). 

The  expected  value  of  J(jr,t)  is  defined  by 

Je(r,t)  =  Tr[J(r)+JA(r,t)]D(t)  (11-52) 

where  D(t)  is  the  density  operator  of  the  perturbed  system  at  time  t  and 


satisfies  the  Liouville  equation 

D(t)  =  ^  [D(t),  HP+H1(t)+H2(t)]  (11-55) 

with  the  boundary  condition  D(t*-oo)=DP. 

£ 

For  a  weak  perturbation,  we  shall  obtain  the  current  response  J/(r,t)  to 
first  order  in  A.  For  this  purpose  we  substitute 

D(t)  =  DP+D1(t)  (11-54) 

into  eq.  ( 11-53 J  where  D^(t)  is  the  perturbation  due  to  A  and  neglect  the  terms 
[Dp,H2(t)]  and  [D1(t) ,H1(t)+H2(t) ]  in  the  resulting  equation  (note  that  H-^  is 
first  order  and  H2  is  second  order  in  A).  Then,  satisfies 

“  DX(t)  =  |  [DX(t),  HP]  [DP,Hx(t)]  (H-55) 

whose  solution  is  readily  found  as 


Bl(t)  =  J 


t 

idt' 


Dp,e 


j:  HP(t-t 1 

.  Tl 


) 


Hi(t')e 


j  HP(t-t 

T1 


]• 


(11-56) 


I 
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The  substitution  of  eqs.  (11-54)  and  (II-56)  into  eq.  ( 11-52)  gives,  after 

2 

neglecting  the  terms  containing  A  (r, t). 


J!(r,t)  =  TrDPJ.(r)  +  TrDPJ,,(r,t) 


+  z  fd^r'  /dt' TrI)P  j  A^'t' 


(H-57) 


where  the  subscripts  £  and  m(i, m=l, 2,3)  refer  to  the  components  of  the  vectors 
and  where  the  summation  convention  on  the  repeated  indices  is  used.  The 
first  term  in  eq.  ( 11-57)  is  the  expected  value  of  the  current  in  the  unper¬ 
turbed  medium  and  is  zero.  The  second  term  is  given  explicitly  by 


™>P4(r,t)  =■  -  E  DS„  <  n  |  6(r-roj)A(roj,t)|n> 


Since  the  delta  function  is  contained  in  the  matrix  element,  the  integration 
over  all  the  coordinates  of  particles  will  give 


p  ,  x  2  £(  r,  ^ ) 

TrD  JA(l,t)  =  -q  - 

2 

where  q  is  defined  in  eq.  ( 11-45).  Then  eq.  ( 11-57)  becomes 


(11-58) 


A  ( r,  t )  +  A  ( r,'  t ' ) 

J^(l,t)  =  -q2  - - -  +  ^  /dVjdt'  TrDP[J£(r,t),  Jm(r,'t')]  —q -  .  ( 11-55 

In  order  to  calculate  the  absorption  and  refractive  index  in  the  frame¬ 
work  of  Kubo's  theory,  we  substitute  eq.  (II-59)  into  eq.  (11-46)  and  use 
eq.  ( 11-47)  to  eliminate  the  vector  potential  in  favo'  of  the  electric  field. 
The  result  is 


2 

(r,t)-  E  (r,t)  =  ~  /d3r’  /dt  ’  {q2P(t-t  ’  )6(r-r'  )E 


c2  dt2  1 
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-  i  TrDP[J  (r,t-f),  J  (r*  )]JE Jrjt' )  (II-60) 

ft  l  m  m 

or 

^E(r,t)-~-  E(r,t)  /d5r'  jdt'  a(r-r,'t-t'  ) -E^t' )  (II-61) 

c2  dt2  c2  -00  =  - 

where  a(R, t)  is  called  the  conductivity  tensor.  We  shall  give  the  explicit 
form  of  a  in  the  transformed  domain  later. 

Equation  (11-60)  which  describes  the  electric  field  in  the  medium  can 
be  solved  if  one  specifies  the  boundary  and  the  initial  conditions  on  E(r,t). 
The  solution  can  in  general  be  constructed  in  terms  of  the  solution  of  the 
homogeneous  equation  of  the  following  form 


3-(ik-r-icut)  ( 11-62 ) 

where  k  and  cu  are  related  to  each  other  through  the  dispersive  relation. 


hnvx> 


.  2  (A)  “Mi  diM  /  v 


where  o  (k,<u)  is  the  scalar  conductivity  defined  by 


(11-63) 


ox(k,u))  =  e^.o(k,o))-€K 

where  a(k,u))  is  the  transform  of  a(R, t)  and  explicitly  given  by 

ox(k,o))  “  "  “J"  -  ^  /d5r'  ^dTe-ia)Tei—  ^TrDP[J^(r, t),  J^(r' )]  (11-64) 

In  the  applications,  one  usually  encounters  two  types  of  problem: 

(l)  initial  value  problem  and  (2)  boundary  value  problem.  In  the  first  case, 

one  specifies  an  initial  spatial  distribution  E(r, o)  and  solves  for  E(r, t) 

for  t>o.  The  initial  distribution  can  be  expressed  as  superposition  of  the 

ik  •  r 

terms  of  spatial  modes  of  the  helmholtz  operator,  i.e.,  e - where  k  is  real 
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vector.  The  damping  associated  with  each  mode  is  obtained  by  solving  the 
dispersion  relation  as  cu(k)  for  a  given  real  k. 

In  the  second  type  of  problem,  one  solves  the  Maxwell's  equation  (II-60) 
for  E(r, t)  when  E(r, t)  is  specified  on  the  boundary  as  a  known  function  of 
time.  The  time  dependence  can  be  expressed  as  the  superposition  of  terms  of 
the  form  e^^  where  o>  is  a  real  number.  The  associated  complex  k  is  obtained 
from  the  dispersion  relation  as  k(cu).  A  sinusoidal  plane  wave  impinging  on 
the  side  of  an  half  infinite  medium  is  a  typical  example  for  problems  of  the 
second  kind.  The  complex  number  k(cu)  in  this  case  is  the  inverse  relaxation 
length  in  the  medium. 

We  shall  compare  the  Kubo  theory  to  transport  theory  in  these  two  typical 
cases . 


a.  Damping  in  Time 

As  we  mentioned  above,  we  must  solve  the  dispersion  relation  in  this 
case  for  a  real  k,  and  obtain  the  real  and  imaginary  parts  of  co(k)  for  a 
given  k.  It  is  convenient  to  substitute 


CD  = 


CD- 


CDr 


(II-65) 


n(k)  n  +in, 

where  cDQ=ck,  and  n(k)  is  called  the  complex  refractive  index  with  nQ  and  n, 
as  its  real  and  imaginary  parts.  The  electric  field  will  decay  in  time  as 
e  where  cdj  is  the  imaginary  part  of  cu(k).  &dj  is  the  decay  rate  of  the 

electric  field  which  is  the  quantity  to  be  compared  to  the  photon  absorption 
coefficient  per  unit  time  a^(k,CD^)  obtained  in  photon  transport  theory. 

The  substitution  of  eq.  (II-65)  into  the  dispersion  relation  eq.  (II-65) 


*  •  , 


2k 


gives 


2  2 
n  -n 


2no"l 


inti 


o.  (k> 


CD 


i-  -2—1-  +  i 

,  2  2  >2  ,  2  2N2  cd  (n  +in  )  V-'  n  +in 

(n0+n1)  (no+nl^  0  o  l'  o  1 


In  order  to  solve  the  dispersion  relation  for  the  real  and  imaginary  parts 
of  ca(k),  we  shall  consider  a  weakly  absorbing  medium,  n  «nQ.  In  this  case  the 
dispersion  relation  can  be  approximated  as 


,  1  ,  2nl  Uni  .  R/,  .  I,,  ,, 

1-  —  +  i  — -  =  -  (k,cDk)  +  i  ov  (k,CDj) 

o  o 


CD 


O  . 


where  cd =  —  is  exactly  the  frequency  used  in  photon  transport  theory  and 
k  nQ 

R  I 

o  (k,cn )  and  a  (k,cD  )  are  respectively  the  real  and  imaginary  parts  of  the 

A.  K  \  K 

scalar  conductivity  a  (k,CD >  )  given  by  eq.  (II-6U). 

\  K 

Equating  the  real  and  imaginary  parts1,  one  obtains  the  real  refractive 

index 

n2  =  1  +  —  a  I(k,cD  )  ,  (11-66) 

o  cd,  \  —  k  ' 

k 

the  imaginary  part 


2nn, 

) 

CDk  \ 


0  X 

n,  =  -  -  (k,cu^) 


and  the  decay  rate  of  the  electric  field 


2cd, 


2nA 


R 


=  Una  (k,cD  ). 


(11-67) 


The  explicit  form  of  a  (k,ca)  was  given  in  eq.  (II-6U)  for  complex  cd. 


CD- 


When  cd  =  —  =  cd,  ,  we  have 
n0  k 


Wtfjr.r),  J^r')] 


By  using  the  integration  representation  of  the  unit  step  function,  i.e.. 
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-iST 


u(t)  =  ir-lds  uith  e"° 


and  letting  t  =  -t'  and  03'  =  s+03 ,  then 

1C 

let  <k,tu  )  -  a!  +  _i_  /dV  7df  /to'  (r),  J  (r}t')]. 

\  k  03  2jta,fi  -So  -So  03-03,+ie  \  -  7  \ 

k  k  k 

The  real  and  imaginary  parts  of  a  (k,03  )  can  be  obtained  by  noting  that 

a.  k 


as 


cvE(S’V 


lim  — — — ■ —  ■  P  —7-  -  in&(tD-ai) 
e-*o  ui-ku+ic  uj^o) 


/dV  /dfe^'^'d-l'Wu  (r;t.),  j  (r)] 

2fto.  -«>  \  \  — 


(11-68) 


*  5^  'd3r'idt'pito'  ^ 

£ 

“k 


00  00 


,,  e 


±03 '  t  ’ 


,ik,(j:-I,)TrDP[J  (p.t.)#  J  (r)] 


_<£ 


(11-69) 


b.  Damping  in  Space 

As  mentioned  above,  the  dispersion  relation  must  be  solved  in  this  case 
for  a  real  0>=o3o  to  obtain  the  real  and  imaginary  parts  of  K(u3q).  It  is  con¬ 
venient  to  let 

K  =  (n  +in. )k 

where  ck=03o.  For  a  weakly  absorbing  medium  n^«nQ,  the  dispersion  relation 
eq.  (II-63)  can  be  approximated  as 


2  4ni  ,  R,  ,  «  ,  1/  , 

n  -1  +  2in,n^  =  -  -  {a,  n  +  io,  (n  k,03  )). 

o  io  03„  \  o->  o'  \  o->  o" 

0 

Then,  the  real  refractive  index  and  the  imaginary  part  are 


(n  k,03  ) 

o  \  0~  ’  O 

n^  =  1  +  - 

0  03_ 


(11-70) 
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and 


2n 

n  cu 
o  o 


\r(^^q)  • 


The  E.M.  wave  absorption  per  unit  length  is 


(H-71) 


\(±Pq)  =  -2nik  “  “T  ^(n^o).  (II-72) 

This  is  the  quantity  to  be  compared  to  the  photon  absorption  coefficient 
per  unit  length  a^k/c^)  obtained  in  photon  transport  theory. 

In  eqs.  (11-70),  (11-71),  and  ( 11-72),  c^R( n£pQ)  and  o^njc,^)  have 
the  exact  forms  as  eqs.  (11-68)  and  ( 11-69)  except  that  the  places  where  k 
and  occupy  are  respectively  replaced  by  n^  and  u>Q. 


5.  Comparison  of  the  Results  From  Two  Theories 

The  resulting  expression  obtained  by  substituting  eq.  ( 11-68)  into  eq.  ( II- 
67)  is  identical  to  eq.  (ll-3l\  and, furthermore,  eq.  (11-66)  with  the  substitution 
of  eq.  (II-69)  is  the  same  as  eq.  (Il-b3)  with  the  substitution  of  eq.  (il-bb). 

One  can  conclude  that  the  photon  absorption  coefficient  per  unit  time  and 

the  refractive  index  obtained  from  photon  transport  theory  are  the  same  as 

that  obtained  from  Kubo's  theory  through  damping  in  time  only  for  weakly  absorbing 

media.  Since  eqs.  (II-36)  and  (11-1*3)  are  not  the  same  as  eqs.  (11-72) 

and  (II-70),  respectively,  the  photon  absorption  coefficient  per  unit  length 

and  the  refractive  index  obtained  in  the  former  theory  are  different  than 

that  obtained  in  the  latter  theory  through  damping  in  space  even  for  a  weakly 

absorbing  medium. 

As  mentioned  in  section  II-l,  the  "dressed  photon"  concept  used  in 
photon  transport  theory  is  to  associate  photons  in  medium  with  the  same 
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2  XT 

wavelength  \  =  —  as  in  free  space,  but  a  different  frequency  from  ck. 

It  is  obvious  that  the  dressed  photon  technique  is  equivalent  to  damping 
E.M.  wave  in  time  because  the  latter  is  to  consider  E.M.  wave  in  medium  as 
a  wave  having  the  same  wavelength  as,  but  a  different  frequency  from  the 
free  space  values.  It  is  different  from  damping  E.M.  wave  in  space  because 
the  latter  is  to  regard  E.M.  wave  in  medium  having  the  same  frequency  as, 
but  a  different  wavelength  from  the  free  space  values. 

From  the  above  conclusions  and  reasons,  it  seems  suggestive  to  for¬ 
mulate  photon  transport  theory  by  associating  photons  in  medium  with  the 
same  frequency  as,  but  a  different  wave  number  vector  from  that  in  free 
space.  Then  the  photon  absorption  coefficient  per  unit  length  and  the 
r  ractive  index  obtained  from  such  formulated  photon  transport  theory 
may  turn  out  the  same  results  as  obtained  from  Kubo's  theory  by  damping 
E.M.  wave  in  space.  It  would  not  give  the  photon  absorption  coefficient 
per  unit  time  and  the  refractive  index  as  obtained  from  Kubo's  theory  by 
damping  E.M.  wave  in  time.  However,  the  suggested  formulation  is  more 
desirable  because  most  measurements  have  been  usually  done  in  measuring 


photon  absorption  coefficient  per  unit  length. 


CHAPTER  III 


VARIOUS  PHOTON  ABSORPTION  MECHANISMS 


In  chapter  II,  we  have  obtained  the  absorption  coefficient  for  a  photon 
with  momentum  fik  and  polarization  \  in  dispersive  media.  We  shall  restrict 
ourselves,  henceforth,  to  a  nondispersive  medium.  In  order  to  display  var¬ 
ious  mechanisms  contributing  to  the  photon  absorption  in  a  material  medium 
and  also  to  estimate  the  order  of  magnitudes  of  the  various  contributions  in 
the  lowest  approximation  we  shall  change  the  representation  {|n>}  for  the 
description  of  the  particle  system.  In  chapter  II,  the  particle  wave  func¬ 
tions  |n>  were  chosen  as  the  eigenfunctions  of  the  total  hamiltonian 


which  included  the  interaction  between  various  constituent  molecules.  In 
the  present  chapter  we  shall  work  with  a  representation  which  diagonalizes 
only  part  of  HP.  The  remaining  part  of  HP  will  be  treated  as  perturbation. 

1.  DESCRIPTION  OF  THE  PARTICLE  SYSTEM  IN  FIRST  ORDER  PERTURBATION  APPROX¬ 
IMATION 

For  a  particle  system  which  consists  of  neutral  atoms,  singly- charged 
ions  and  free  electrons,  the  hamiltonian  of  the  system  can  be  written  as 

HP  =  ^  +  H1  +  He  +  VP 

vp  .  +  vee  +  v11  +  /e  +  Z1  +  vie 

where  Z,  Z,  and  H6  are,  respectively,  the  hamiltonians  of  the  neutral 
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atoms,  the  singly- charged  ions  and  the  free  electrons.  V  ,  V  and  Vx 
are  the  potentials  between  the  same  kind  of  "molecules"  (atoms,  ions  and 
free  electrons),  and  Ze,  Z1  ,  and  V*e  are  the  potentials  between  the  dif¬ 
ferent  kind  of  "molecules." 

In  addition  to  H^,  the  interaction  of  the  particles  with  the  radiation 
field  can  be  separated  as 

HpR1  _  v  +  v  +  v 

A  e  i 

where  V^,  Vp,  and  are,  respectively,  the  interactions  of  the  atoms,  the 
free  electrons,  and  the  ions  with  the  radiation  field. 

A  convenient  and  simple  representation  which  one  can  choose  is  that 
the  wave  functions  |n>  =  |aP7>  satisfy  separately  the  SchrBdinger  equations 

HA|a>  =  Ea|a> 

He|P>  =  Eg|p> 

h1  | r  >  =  E7|7> 

p  pRi 

with  V=V  +  H  considered  as  the  perturbed  hamiltonian.  Since  each  of 
h\  H  ,  and  can  be  separated  as  the  sum  of  the  hamiltonians  for  the 
molecules  of  the  same  kind,  the  wave  function  |aS7>  is  the  product  of  the 
wave  functions  for  the  individual  molecules  in  the  system.  In  this  rep¬ 
resentation  we  must  calculate  the  transition  probabilities  at  least  to  the 
second  order  in  V  if  we  want  to  investigate  absorption  due  to  free-free 


transitions  of  charged  particles.  The  photon  transport  equation  then  becomes 
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^fxfck,t)+  £'Wf^(r,k,t)  *  I! — (>i-T|)|<n,T1,|v|ni|> 


+  V  <n'i'|v|i."n"><n'V|T|mi>  ,g  ■  , 

„V  E  -E  ,  ,  1  MErf,fV 


(m-i) 


For  the  reduction  of  eq.  (ill-l),  we  shall  use  second  quantization 
in  which  the  potentials  /*,  Z1,  and  Vie  and  the  interactions  VA>  V  t 


and  are  given  by  (see  Appendix  B) 


2 

Ap  4ne  &k(K-K'+u-u '  )  2  w.-.-tx. 

v  =  I - - - - - <a'|-z+Xe - ^j|a>At(K'a' )At(u' )A(Ka)A(u) 

M'  L3|u-u'|2  Jal 


(III-2 ) 


Z1-  Z  ,z_  l  e-i(K-ie'  )-fj  ,,  ><b.  |z.Vel(3S-K'  )~ft  ,b> 


1 l5|k-k'|2 


a  A^(K'a'  )A+(_i,b'  )A(Ka )A(jb) 


(HI-3) 


Une  &k(jL-1'+u-u'  ) 


-z  +Zi1ei^-"-'  )^j|b>At(/'b'  )A+(u')A(_ib)A(u) 

J 


(III-U) 


VA  =  &K(K-K'-is)^a.a £a-a  •[c^+(k)€.K(k)+c^(-k)<€x(-k)]At(K'a’  )A(Ka) 

Ka  L  0) 

]Ta'  (III-5) 


2 

ve  *  2*T  6K(u-u'-k)A+(u' )A(u) 

-  ,  h(X) 


(III-6) 


1 


r-t  1  G 

v.  =  Z  — 

1  Xk  c 

ib 

i'b' 


where 


gn-hc 

t3 

L  u> 


—  6K(i-;,-k)cub!bdb,b  -{C^(k)€x(k)+C^(-k)ex(-k)}  A^/’b'jAQb) 

(HI-7) 


da, a  =  <a'l  ZrP(j|a>  ,  dfe,b  =  <'b’|^£j|b>  ,  <V 


E  ,-E 
x  x 
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In  the  above  equations,  A(u),  A(Ka),  and  A(/b)  are,  respectively,  the  destruction 
operators  for  the  annihilation  of  a  free  electron  of  the  momentum  -ftu,  an  atom 
of  the  external  (center-of-mass )  momentum -fiK  and  the  internal  state  |a> 
and  an  ion  of  the  external  momentum  fi/  and  the  internal  state  |b>.  A^(u'), 

A^(K'a'),  and  A*(_/'b')  are  the  creation  operators  which  create  a  free  electron, 
an  atom  and  an  ion  in  different  external  and  internal  states.  The  symbol  & 
with  subscript  K  denotes  the  Kronecker  delta.  is  the  position  of  the  j-th 
atomic  electron  in  an  atom  or  ion  with  respect  to  the  position  of  its  nucleus. 
e^a'a  (°r  e^b'b^  the  rrabr^x  element  of  the  dipole  operator  with  respect 
to  the  internal  atomic  (or  ion)  states.  '^)a>a<a  (°r  "^b'b^  bbe  energy 
difference  of  the  internal  states  |a'>and  |  a  >  (  or|b' >  and|b^  of  an  atom 
(or  ion). 


2.  TRANSITION  PROBABILITY  FOR  PHOTON  EMISSION 

In  this  section,  we  shall  make  several  approximations  to  reduce  eq. 
(ill-l)  to  a  simple  form.  The  number  of  photons  with  momentum  dik  and  po¬ 
larization  X  emitted  from  l?  per  unit  time  can  be  obtained  by  letting  rj ' = rj+ 1 
in  eq.  (III-l),  i.e. 
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E  (k)(f  (k)+l)  =  Z^D  <n,Ti+l|v|nr1>+  Z 

^  nn,  *  nr), mil  i'll  n%" 

q  =nq,n'q+l 


<n '  q+1  V  n  n  >  <n  n  V  n 


E  -E  it  it 
nrj  n  q 


®  5(E  i  . ,  -E  )  (III-8) 

n  r  l  nq 

where  E^(k)  is  the  transition  probability  per  unit  time  for  emission  of  a 
photon  of  momentum -fik  and  polarization  Let  us  consider  first  the  direct 
transition,  the  first  term  in  the  absolute  square  of  eq.  (III-8).  Since  the 
potential  V*5  between  the  particles  contains  neither  the  photon  creation 
nor  destruction  operator,  then 


<n'T]+l|  V|nr)>  =  <a'  q+l|vA|aq>  5pp  <by7<  +  '  q+l| Ve | Piq >  6^'  &77  • 

+  <7  '  q+llvj  7q>  i 

where  <p '  q  +l|ve|Pq>  '  accounts  for  the  direct  photon  emission  through 

the  interaction  of  the  free  electrons  with  the  radiation  field.  No  such 
transition  can  exist  because  the  energy  and  momentum  conservation  laws  can 
not  simultaneously  hold  (see  the  justification  in  Appendix  C).  Then 


<n'  q+l|  V|  nr)  >  =  <&'  q+l|VA|aq>  ,6^,  +  <7  '  T+l  I  \  1 7h>  5^  .  . 

The  second  term  in  the  absolute  square  of  eq.  (III-8)  accounts  for  the 
indirect  photon  emission  through  the  intermediate  states  |n"q">  .  It  has 
the  non-zero  contribution  only  for  q"=q  and  q"=q+l. 


n"t|"/nr),n'  t)+1 


<n'n+l|vln"n"  ><n"n' 

E  -E  „  ,, 
nq  n  q 


lvlnn>  _  -  <n '  I VP I  n"  >  <tl"n+l  I 

W 


33 


where  we  have  used  E  ,  ,-E  =o  from  the  energy  conservation  delta  function 

n'Tj+l  np 

in  obtaining  the  demoninator  of  the  first  sum.  To  simplify  eq.  (ill- 10), 
the  following  three  assumptions  will  be  made  throughout  this  chapter. 

(i)  The  neglect  of  the  potentials  between  the  molecules  of  the  same 
kind  (i.e.,  neglect  of  V(  e,  and  V11).  This  assumption  is  justified 

because  these  potentials  do  not  affect  the  photon  emission  and  absorption 
much,  although  they  play  an  important  role  in  the  shape  of  the  lines  (pres¬ 
sure  broadening).  Under  this  assumption,  eq.  (III-10)  becomes 
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The  first  sum  accounts  for  the  photon  emission  due  to  electrons  moving 
in  the  field  of  neutral  atoms.  We  shall  refer  to  this  as  the  bremsstrahlung 
of  electrons  in  the  field  of  atoms.  This  sum  consists  of  two  parts.  The 
first  one  corresponds  to  the  emission  mechanism  in  which  an  electron  interacts 
first  with  an  atom  through  the  coulomb  potential  and  goes  into  an  intermediate 
state.  It  then  interacts  with  the  radiation  field  by  emitting  a  photon. 

The  second  part  represents  the  emission  mechanism  in  which  an  electron  first 
emits  a  photon  and  then  interacts  with  an  atom.  The  second  sum  in  eq. 

(III-11)  is  the  bremsstrahlung  of  electrons  in  the  field  of  ions  and  the 
third  sum  is  the  bremsstrahlung  of  ions  in  the  field  of  atoms.  Each  of  these 
sums  contains  two  parts  corresponding  to  the  same  emission  mechanisms  as 
described  above. 

The  fourth  (fifth  and  sixth)  sum  is  the  dipole  radiation  of  atoms 
(atoms  and  ions)  induced  by  electrons  (ions  and  electrons).  The  reason  we 
use  this  terminology  is  that  the  emitted  radiation  comes  from  the  dipole 
transition  of  the  atoms  induced  by  the  interacting  electrons. 

(ii)  The  internal  states  of  the  ions  are  unchanged.  This  implies  that 
V^=o  because  u>^i^=o  for  b'=b  (see  eq.  (III-7)).  Then  the  third  and  sixth 
sums  in  eq.  (ill-ll)  vanish  and  eq.  (III-8)  becomes 


<n ' tj+  1 1 V | rnq >  =  <a' q+l|  VA|aq>  5pp  1 6^' 


(HI- 12) 
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under  this  assumption. 

For  structureless  ions,  such  as  hydrogen  ions,  Vj  takes  the  form  of 
eq.  (III-6)  for  Ve.  In  this  case,  the  third  and  sixth  sums  can  be  neglected 
as  compared  to  the  first  two  sums  because  the  mass  of  hydrogen  ion  is  much 
larger  than  the  electron  mass.  Furthermore,  eq.  (III-12)  is  also  true  for  this 
case  because  the  free  hydrogen  ion,  just  as  the  free  electron,  can  not  emit 
or  absorb  any  photon  through  direct  transition. 

(iii)  All  the  cross  terms  after  expanding  the  absolute  square  of  eq. 
(III-8)  will  be  neglected. 

Under  the  above  three  assumptions,  eq.  (III-8)  becomes 
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are  the  numbers  of  photons  of  momentum  hk  and  polarization  \  emitted  per  unit 

3  BA  • 

time  from  L  ,  respectively,  due  to  the  atomic  bound-bound  transition  (e  ), 

A. 

/  BrAe \ 

the  bremsstrahlung  of  electrons  m  the  field  of  atoms  (e  )  and  ions 

A. 

(e  Brie)^  the  atomic  dipole  moment  transitions  induced  by  electrons  (e  °^e) 

A.  A. 

and  ions  (c  El^)* 

Since  the  radiation  absorption  due  to  electrons  in  the  field  of  atoms 
will  be  investigated  in  Chapter  IV,  the  reduction  of  eq.  (III-13)  to  the 
form  used  in  Chapter  IV  is  now  performed  in  detail. 

Let  the  initial  and  final  states  of  the  electrons,  the  atoms  and  the 


photons  in  the  system  be 


|af3->  =  |...n(Ka),  n(K' a |  ,n(u),  n(u'),  n(u  ")...> 
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|q+l>  ■  |  . .  .T^(k)+l. . .  > 

where  n(Ka)  is  the  occupation  number  of  the  atoms  of  the  external  momentum  <fik 
and  the  internal  state  |  a  >.  n(u)  is  the  occupation  number  of  the  electrons 
of  momentum  <fiu  and  same  for  T^(k).  One  the  intermediate  states  of  the 
electrons  for  which  the  matrix  elements  in  eq.  (III-15)  does  not  vanish  is 


|p">  =  |...n(u)-l,  n(u'),  n(u")+l...> 

From  eqs.  (III-2)  and  (III-6),  eq.  (III-I5)  becomes 
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N(u),  N(Ka)  and  f  (k)  are  the  numbers  of  the  electrons,  the  atoms  and  the 
photons  in  l? ,  respectively.  In  obtaining  the  above  equation,  n(u)  and 
n(Ka)  have  been  neglected  in  comparison  with  unity  and  the  statistical  average 
over  the  initial  particle  and  photon  states  have  been  performed.  For  non- 
relativistic  electrons  one  can  replace  Bjj(u-ulk)  by  6k(u-u').  It  means 
that  the  recoil  momentum  of  the  electrons  can  be  neglected.  Taking  the 
sum  over  u"  after  this  approximation  is  made,  one  obtains 
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where 


h  5 ( K-K '  +u-u  '  )&  ( E  ,  ,  -E  4-E  ,-E  4*u)  (III-22) 

i\,  &  1Y&  U.  U 

=  u-u'  and  =  £»e  . 

N.  and  N  are  the  number  densities  of  the  atoms  and  electrons  in  the  system. 

A  0 

P  =N(a)|N  is  the  ratio  of  the  atomic  density  in  the  internal  state  |a>  to 
the  total  density  of  the  atoms  in  the  system.  It  can  be  interpreted  as  the 


39 


probability  of  finding  the  atoms  in  the  state  |a>.  M^(K)  and  Mp(u)  are 

3  3 

the  distribution  frequencies  of  the  atoms  in  d  K  and  of  the  electrons  in  d  u. 
The  further  reduction  of  eq.  (III-22)  will  be  found  in  Chapter  TV. 

In  the  same  way,  the  transition  probabilites  per  unit  time  for  photon 
emission  due  to  the  other  mechanisms  as  described  in  eqs.  (ill-lt),  (lV-16), 
(III-17),  and  (III-18)  can  be  written  down.  Letting  the  ion  density  of  the 
system  be  equal  to  the  electron  density  the  total  transition  probability 
per  unit  time  for  emission  of  a  photon  of  momentum  fik  and  polarization  \ 
due  to  all  the  significant  mechanisms  under  the  assumptions  made  is 
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l^'a'41  *  Vb(3>  -  i<b'i-J^2-J|b>l  • 

In  eqs.  ( II 1-25 )  and  (III-27),  is  the  probability  of  finding  an  ion 

in  the  state  |b>  and  describes  the  distribution  of  the  ion  velocities. 


3 .  PHOTON  ABSORPTION  COEFFICIENT 

It  has  been  shown  in  section  III-2  that  the  transition  probability 
per  unit  time  for  emission  of  a  photon  of  momentum  *k  and  polarization  X 
due  to  all  the  significant  mechanisms  E  (k)  is  given  by  eq.  (III-23)-  It 

K  ■*“ 

is  possible  from  E  (k)  to  calculate  the  absorption  coefficient  for  photons. 
K 

This  calculation  is  now  performed. 

If  '.i.  is  assumed  that  the  medium  is  isotropic,  then  the  absorption 
coefficient  for  unpolarized  photons  is 


(III-28) 


where  A  (k)  is  the  transition  probability  per  unit  time  for  absorption  of 

K 

a  photon  of  momentum  fik  and  polarization  X. 

For  M^(K),  (u)  and  being  Maxwellian  distributions  and  for  the 

internal  states  of  the  atoms  (ions)  populated  as  e  ^  (e  where  6  is 


iil 


is  the  temperature  of  the  system  in  energy  unit,  then  the  absorption  and 

emission  transition  probability  are  related  by  the  equation 

fix 

\(k)  =  e  6  Ex(k).  (III-29) 

From  eqs .  (111-23),  (III-28)  and  (III-29),  the  absorption  coefficient  can 
be  written  as 


4.  TIME-DEPENDENCE  OF  a 

A  plasma  will  vary  with  time  and  eventually  die  out  if  there  are  no 
external  devices  to  maintain  it.  The  absorption  coefficient  a  is  a  function 
of  time.  The  possible  parameters  in  eq.  (III-30)  which  may  depend  on  time 
are  the  neutral  and  electron  densities,  the  temperature  of  the  medium  and 
the  probabilities  Pa  and  P^  of  finding  respectively  the  atoms  and  the  ions 
in  the  states  |a>  and  |b>.  Since  plasma  temperatures  are  known  to  be 
sensitive  functions  of  time  for  most  cases,  CQ,  and  are  time-dependent. 


h2 


In  this  section,  we  shall  investigate  the  time  variations  of  the  quantity 
y(t)  =  - instead  of  a(t),  by  assuming  that  the. neutral  and  electron 

NeoCi(t) 

densities  satisfy,  respectively,  the  simple  differential  equations 
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70Ni(t) 


(III-3M 


=  ,i(,) 

dt  0  e 


(III-35) 


where  yQ  is  the  recombination  coefficient  of  the  electrons  with  the  ions. 
The  solutions  of  eqs.  (III-3U)  and  (III-35)  are 


Neo 
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where  Ngo  and  N^Q  are  the  respectively  the  neutral  and  electron  densities  at 
the  instant  of  the  plasma  formation.  The  substitution  of  eqs.  (III-36)  and 
(III-37)  into  eq.  (III-30)  gives 


y  =  +  &+&>  (m-38) 

where  x  =  r0Neot,  y(t)  =  a(t)/N^QC1(t) ,  \  =  NAo|Neo  and  |i(t)  =  C2(")/C1(t) 
are  dimensionless  positive  numbers.  In  obtaining  eq.  (lIx-38),  the  di¬ 
rect  bound-bound  transitions  of  the  atomic  electrons,  i.e.,  the  first  term 
in  eq.  (III-30)  has  been  neglected.  The  reason  for  neglecting  this  term 
is  that  it  represents  the  atomic  line  absorption,  and  hence  is  negligi¬ 
ble  when  the  frequency  of  the  photons  is  far  away  from  the  line  fre¬ 
quencies.  In  eq.  (III-38)  \  is  the  ratio  of  the  neutral  to  the  electron 


^3 


density,  a  quantity  of  measuring  the  degree  of  ionization  of  the  medium  at 
the  instant  of  the  plasma  formation,  n  is  the  ratio  of  the  inverse  brems- 
strahlung  due  to  an  electron  in  the  field  of  an  ion  to  that  of  an  electron 
in  the  field  of  an  atom  if  the  contributions  due  to  induced  dipole  transitions 
are  negligible  (see  eq.  (III-32)).  In  the  case  that  C^(t)  and  C?(t)  are 
sensitive  to  the  time  t,  their  ratio,  p(t),  may  be  insensitive  to  t  because 
the  numerator  and  the  denominator  are  both  time  dependent  through  the  tem¬ 
perature  in  the  Maxwellian  distribution  of  particles.  We  shall  assume  that 
H  is  constant,  then  a  maximum  of  y  in  eq.  (111-38)  occurs  at  a  value 
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provided  \+2p<l  because  \,  n  and  x  are  positive  numbers.  The  variation 
2 

of  y(t)=a(t  )/neQC-L(t)  with  time  is  shown  in  Figure  1  for  the  plasmas  with 
\=l/9  and  p=o,  l/9,  3/9,  3/9  and  b/9.  If  the  dependence  of  plasma  tem¬ 


perature  with  time  is  known,  the  variation  of  the  absorption  coefficient 
a(t)  can  be  obtained  through  the  calculation  of  C^(t). 


CHAPTER  IV 


RADIATION  ABSORPTION  IN  PARTIALLY  IONIZED  HYDROGEN 
GAS  DUE  TO  INVERSE  BREMSSTRAHLUNG  OF  ELECTRONS 
IN  NEUTRAL  ATOMS 

(12) 

In  1967  Akcasu  and  Wald  investigated  the  radiation  absorption  due  to 
the  inverse  bremsstrahlung  of  slow  electrons  in  the  field  of  neutral  atoms. 
Since  the  temperature  of  the  system  they  investigated  was  low  (©-lev  or  less), 
they  assumed  that  all  the  atoms  in  the  system  are  in  the  ground  state,  the 
energies  of  the  electrons  are  insufficient  to  excite  an  atom  from  its  ground 
state  to  an  excited  state,  and  the  elastic  scattering  cross  section  for 
electron-atom  collisions  appearing  in  the  absorption  formula  can  be  approx¬ 
imated  by  its  value  at  zero  electron  energy.  Under  these  assumptions  they 
calculated  the  various  absorption  contributions  due  to  inverse  bremsstrahlung, 
induced  dipole  transition,  and  exchange  and  interference  effects;  and  found 
that  the  last  three  contributions  for  low  temperature  system  are  negligible 
as  compared  to  the  first  one. 

( 13  ) 

For  hot  plasmas,  such  as  the  one  produced  by  a  giant  pulsed  laser 
beam  which  we  shall  discuss  later,  the  temperature  of  electrons  in  the  plasma 
is  about  20  eV.  At  such  temperatures,  the  above  assumptions  made  by  Akcasu 
and  Wald  cannot  hold.  It  is  the  aim  of  this  chapter  to  consider  the  problem 
for  higher  electron  energies.  The  atoms  in  the  system  are  allowed  initially 
and  finally  to  be  in  any  excited  state  as  we]l  as  in  the  ground  state.  The 
electron  energy  dependence  of  the  elastic  and  inelastic  cross  sections  will 
be  also  taken  into  account.  However,  only  the  absorption  due  to 
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the  inverse  bremsstrahlung  of  high  energy  electrons  in  the  field  of  excited 
and  ground  atoms  shall  be  computed.  The  other  contributions,  such  as 
induced  dipole  transition,  etc.,  which  at  high  temperature  might  not  be 
small  as  predicted  at  low  temperature,  will  not  be  considered  in  this  thesis. 

1.  ENERGY  INTENSITY  OF  EMITTED  RADIATION 

In  section  III-3,  we  have  obtained  in  eqs.  (III-21)  and  (III-22 )  the 
transition  probability  per  unit  time  for  emission  of  a  photon  of  momentum 
fik  and  polarization  \  due  to  the  bremsstrahlung  of  the  electrons  in  the 
field  of  atoms  as 

,  Qa'a(q)  2 

5u'  I  N(a)  Z - 7 -  q  M  (K) 

a  a '  A  - 

■  B(K-K'+u-u' )8(E  ,-E  +E  ,-E+ftu)  (iV-l) 

K  A.  K.A  U  U. 

where  N(a)=N  P  is  the  number  of  the  atoms  in  the  internal  state  |a>  .  In 
A  a 

eq.  (IV-l)  we  have  dropped  out,  for  the  time  being,  the  integration  Ne/d\Me(u) 
which  accounts  for  the  effect  of  the  Maxwellian  electron  distribution  in  order 
to  simplify  the  writing  of  the  expressions  below.  We  shall  resume  this 
integration  later  in  section  (IV-U).  The  superscript  BrAe  in  eq.  (III-21) 
for  indicating  the  contribution  due  to  bremsstrahlung  of  the  electrons  in 
neutral  atoms  is  also  dropped  out  for  the  same  reason. 

For  an  isotropic  medium,  the  energy  intensity  per  unit  energy  emitted 
in  all  directions  and  in  two  polarizations  is  related  to  E  (k)  by  the  equa- 


\(k)  =  /d5 K'  fd 


( IV- 2) 


L7 


S(ftu)  = 


a) 


,  T 

c  (2nr  *• 


where  is  the  direction  in  which  the  photon  is  emitted.  Recalling  that 
q  =£•£  ,  one  has 

A.  A. 


t,  ,  y  .  8  ft  2 

1%  l  \  -  T  q  ’ 


and  eq.  (IV-2)  becomes 


S(-hu)  =  ~ ~~  /  dK5/d5K'  Jd5u*  q2  £  N(a)  F  ,  (q)M  (K)B(K-K'+u-u  ’  ) 

3nmc3  80 '  9  8 


6(Ek-Ek,+A) 


where 


and 


Fa.a(s)  =  l<a' I z- t  e1£'-^|a  >  | ^/q4 

J  — 


(IV-3) 


A  =  E  -E  ,+E  -E  ,-fuj. 
a  a  u  u 

For  a  medium  in  thermal  equilibrium,  the  atom  momenta,  fik,  are  distrib¬ 


uted  according  to  the  Maxwellian  distribution  law,  i.e., 


maU)-j'k  .  -J~ 

(2iiM0r/2 


2  2 

ft  K 


e"  dR T  d  K 


(IV-1, 


3  3 

and  the  integrations  over  d  K  and  d  K1  in  eq.  (IV-3)  can  be  carried  out 


8e 


6  _  _  r?  X  .A  .  m  2/ 

S  (fta>)  =  — - -  /dV  q  L  N(a)F  ,  (3)  p’  ^  Fu  W  Y  ' 

3nm  <T  aa’  8  8  -ftqv2n0 


(IV-5) 


ME 

u  q  3  m 

where  x  =  — 77  and  y  =  .  Since  d  u ’  =  -7s  u'dE  -dQ  ,  ,  the  integration  over 

mQ  u  ft-  u  — u 

cD  ,  can  be  performed  and  eq.  (IV-5)  reuuces  to 


1  p 

,  —  x  ,  A  m  2.c 

S(fto)  »  —  a  -  /dE  ,  /q™ax  dqq2  1  N(a)F  ,  (q)(-)?  e"  Uy?  'Eu  '  M  y  ’ 

5  -2  u  Vn  »»'  8  8  ” 


fiu 


(IV-6) 


where  a  =  -z~  is  the  fine  structure  constant,  q  .  =  lu-u'l  and  q  =  u+u' 

/p  '  Yt  1  1  rttOV 


min 


max 


M 


As  a  result  of  the  large  mass  ratio  -  of  an  atom  to  an  electron,  the 

m 

quantity  x  is  a  large  number  in  most  cases.  Therefore  one  can  approximate 
the  exponential  factor  in  eq.  (IV-6), 


/ x; p  ,  x  .  A  m  2.^.  ,  „  /  A  my 

<;>  exp(-  ~2  %  -  s y  >  >- by  *>• 


The  last  term,  my/  2M,  in  the  argument  of  this  delta  function  accounts  for  the 
recoil  energy  of  th<_  atom.  Thus  eq.  (IV-6)  includes  the  effect  of  finite 
atom  mass.  However,  for  the  sake  of  simplicity,  the  atoms  in  the  medium  will 
be  assumed  to  be  infinitely  heavy  for  the  remainder  of  this  chapter.  Then 
eq.  (IV-6)  becomes 


3  mu 


JdEu'  /qmaxdqq5  EN(a)F  (q)6(E-E  +E  -E  -to) 
o  ^min  aa'  a  a  a  a  u  u 


( IV- 7 ) 


2.  DIFFERENTIAL  CROSS  SECTION 

In  this  section,  we  shall  express  the  energy  intensity  of  emitted 
radiation  in  terms  of  the  differential  cross  section  for  electron-atom 
collision.  During  collision,  the  electron  momentum  changes  from  -ftu  to-ftu' 
while  the  atom  simultaneously  undergoes  a  transition  between  the  initial  and 
final  states  |a>  and  |a’>  for  which  the  atom  internal  energies  are  E  and 

Q 

E  , ,  respectively.  The  differential  cross  section  in  the  laboratory  system 

( liTf  18 ) 

of  coordinates  for  such  a  process  can  be  defined''  *  as 
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a  ,  (E  ,E  ,,q)  dfi  ,dE  , 
aa'  uJ  u  ’  -u'  u' 


O  N  T7  1  o  »  * 


(IV-8) 


(q)6(E  -E  ,+E  -E  .Jdfl.dE 


a2  Eu  a  a  '  a  a'  u  u '  ’  — u '  u 


where  dfi  ,  is  the  element  of  the  solid  angle  in  the  direction  of  u'  and 

— u  - 

■fi2 

a  =  -  is  the  first  Bohr  radius  of  hydrogen  atom.  The  macroscopic  dif for¬ 

me2 

ential  cr  's  section  for  the  scattering  process  can  be  defined  as 


Z  ,  (E  ,E  , ,q)  =  N(a)o  ,(E  ,E  , ,q) 
aa'  u*  u'*  aa'  u'  u'* 


(IV-9 ) 


Since  the  various  atomic  transitions  accompanying  scattering  processes 
remain  unseparated  experimentally,  the  macroscopic  differential  cross  section 
is  obtained  by  summing  the  contributi  vr.s  of  ti  ansi  ‘dons  to  all  admissible 
final  states  of  the  scattering  atom  as 


Z  (E  ,E  , ,q)  =  ~'/rr-  Z,  N(a)F  .  (q)b(E  -E  ,+E  -E  ,  ) 
a  uJ  u  ’  a  2  Eu  a'  a  a  a  a  u  u' 


Then  the  total  macroscopic  differential  cross  section  of  an  electron  scattered 
by  the  atoms  in  all  possible  initial  states  is 


u'  u‘ 


U  Z  N(i 


Z  (E  ,E  ,q)  =  -'J-T-  Z  ,  N(a)F  .  (q)6(E  -E  ,+E  -E  .  ). 
-i'11'  2  Eu  aa’  a  a  «  »  ”  » 


a  a  u  u 


(IV- lx) 


In  terms  of  the  differential  cross  section,  eq.  (IV-7)  becomes 


S(fto)  =  Z  S  ,(ftx) 
aa’  aa 


(IV- 12a) 


wliere 


Saa'(to) 


4  j  2 


a" a  —  jdE  ,  •J—  ■  ~  dqq^  Z  ,(E  ,E  ,  +’ftu,q) 


o  mu  6  u’  Eu,4fu)'qmin 


aa  u  u 


( IV- 12b ) 


i  t 


* 


S(ftD)  «  i  0?a2  ^  /dE  ,  XrV  /qmaW  E  (E  ,E  4*o,q). 
3  o  mu  o  u'  E  ,-tfto)  q  .  u*  u  >H/ 

u  u 1  Tnin 


(IV- 12c) 


If  all  the  atoms  in  the  medium  are  initially  and  finally  in  the  ground 
state  as  assumed  by  Akcasu  and  Wald,  eqs.  (IV-7)  and  (IV-8)  become  in  this 


SM  “  NA  3  ^  mu  ^Eu’  •fqmindqq?F(q)&(Eu-Eul-hcu)  (lV-13) 


o(q)  =  — “  F(q) 

ao 


(IV- 14) 


where  F(q)  denotes  the  matrix  element  evaluated  with  respect  to  the  ground 
state  of  the  atom  and  is  the  number  density  of  the  atoms  in  the  system. 
Equation  (lV-14)  is  just  the  microscopic  differential  cross  section  of  an 
electron  elastically  scattered  by  an  atom  in  the  ground  state.  In  terms 
of  the  velocity,  v=ftu|m,  of  the  electron,  the  integral  microscopic  cross 
section  can  be  written  as 


~  I1  F(q)d|i  =  M“”)2  /  *  qF(q)dq. 
a2  -1  a0mv  0 


(IV- 15) 


Combining  eqs.  (lV-13)  and  (IV-15)  gives 


„  16  5  2  ,m,5  1  ,°°  ,vmaxJ  3  r  /  \  v  da  - 

s(to)  -  N.  —  orV  (-)  -  /dE  ,  /  dw  [ff(v)  +  -  —  j 


A  3tc  o  'fi'  u  v0  u'  J  v 


b  &(E  -E  ,-f/ii) 
u  u' 


(IV- 16) 


V+V 1  I y-v *  I 

where  vmnv  =  “ —  and  vm«n  =  •  Eq.  (lV-l6)  indicates  that  the  intensity 


of  the  bremsstrahlung  is  not  determined  by  the  value  of  elastic  scattering 
cross  section  at  the  incident  electron  energy,  as  might  be  expected  in¬ 
tuitively,  but  it  depends  on  the  variation  of  the  cross  section  in  the 


velocity  region  (v  ,  -v  ).  If  a(v)  is  slowly  varying  up  to  the  incident 
min  m&x 

energy  of  electrons,  one  can  approximately  evaluate  the  integral  over  dv 


in  eq.  (IV-6)  for  a  low  temperature  system 


\  Mi  /  \  /  U\3/2,  -ftuK ,,  ilu).2 
S(hiu)  =  Na—  a(o)c(— )  '  (2-  —  )(l-  a 

me  u  u 


(IV- 17) 


where  a(o)  is  the  cross  section  of  the  electron  elastically  scattered  by  the 


atom  in  the  limit  v-*o.  Equation  (IV- 17)  is  identical  to  that  obtained  by 

(12) 

Akcasu  and  Wald  '  by  the  partial  wave  method.  This  identity  shows  that  the 
method  of  this  chapter  by  using  plane  wave  for  the  electron  wave  function 
will  yield  in  the  above  approximation  the  same  result  as  obtained  by  partial 


wave  method  if  one  uses  the  experimentally  measured  scattering  cross-section 
of  ground  state  atoms  in  both  methods . 

Although  eq.  (lV-l6)  is  obtained  by  assuming  that  all  the  atoms  in  the 
medium  are  initially  and  finally  in  the  ground  state,  it  is  also  applicable 
to  the  atoms  being  initially  and  finally  in  the  same-  excited  state.  In 
this  case,  a(v )  in  eq.  (lV-l6)  is  the  cross  section  of  an  electron  elas¬ 
tically  scattered  by  the  atoms  in  the  excited  state  and  N  is  replaced  by 

A 

N(a),  the  number  density  of  the  atoms  in  the  excited  state  |a>  . 

For  inelastic  scatterings  of  electron-atom  collision,  the  energy  in¬ 
tensity  emitted  through  such  processes  can  also  be  obtained  from  eqs. 

(IV-12)  by  knowing  the  differential  inelastic  scattering  cross  sections. 
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_ 

It  is  obvious  that  the  energy  intensity  calculated  from  eqs.  (IV-12) 
and  (IV-16)  will  be  more  accurate  if  the  experimentally  determined  cross 
sections  are  available.  Unfortunately,  the  cross  sections  are  not  all 
experimentally  measured.  We  shall  calculate,  for  the  sake  of  consistency, 
all  the  relevant  cross  sections  in  Born  approximation  even  though  some  of 
them  have  been  experimentally  determined.  In  order  to  get  various  cross 
sections  in  Born  approximation,  one  has  to  calculate  Faia(q)  for  various 
atomic  states.  F&ia(q)  is  also  contained  in  eq.  (lV-7)>  the  expression 
for  the  energy  intensity  of  emitted  radiation.  However,  once  F  (q)  is 
calculated,  one  can  obtain  the  intensity  of  the  emitted  radiation  directly 
using  eq.  (IV-7)  rather  than  first  evaluating  the  cross  section  and  then 
using  eqs.  (IV-12)  and  (IV-16). 

As  we  shall  apply  the  theory  developed  above  to  the  hydrogen  case  in  the 
next  section,  it  is  suitable  here  to  show  some  of  the  electron-hydrogen 
cross  sections  calculated  in  Born  approximation  and  partial  wave  method. 

The  elastic  cross  sections  calculated  in  Born  approximation  when  the  hydrogen 
atoms  are  initially  and  finally  in  the  ground  state  |l00>  as  well  as  in  the 
states  |200>,  |  210  >  ,  |  211  >  ,  and  |  21,-1  >  of  the  first  excited  energy 
level  are  given  by 
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a  a 

211  21- 


i.-S- 

2  2 
C  RI!  V_ 


l*  3 
10b  +50b 


rran  5U 


b^+lOOb  +95b+37~] 

(n-b)5  J 


2  2  2  2  2  2 

where  b=lU  and  U  =u  aQ  with  -h  u  /2m  being  the  incident  energy  of  the 
electron.  These  cross  sections  are  plotted  in  Figure  2  together  with  a* 
which  is  calculated  by  partial  wave  method ^  and  agrees  well  with 


experimentally  measured  results 


Oioo  was  calculated  before  by  Mott 


and  Massey'  ,  but  for  the  states  in  excited  levels  no  elastic  cross  sections 


for  electron-hydrogen  collision  have  been  calculated  in  either  Bom  approx¬ 


imation  or  partial  wave  method.  Therefore,  we  had  to  calculate  agoo*  ^210 
and  02H  using  the  Born  approximation  explicitly  for  the  purpose  of  compar¬ 
ison  with  o10q.  For  inelastic  cross  sections  when  the  hydrogen  atoms  are 

initially  in  the  ground  or  in  an  excited  state,  most  calculations  in  the  lit- 

( 19  22-27) 

erature  have  been  performed  on  Bc.n  approximation.  7>  ' 

Figure  2  shows  that  the  elastic  scattering  cross  section  for  the  hydrogen 

p 

atom  in  the  ground  state  cal-ulated  in  Born  approximation  is  less  than  o-^qq 
which  agrees  well  with  the  measured  values.  The  discrepancy  becomes  large  when 

the  incident  energy  of  electron  decreases.  Furthermore,  the  elastic  cross  sec¬ 
tion  for  the  atom  in  a  higher  level  is  much  larger  than  that  for  the  atom  in 

a  lower  level  because  the  size  of  the  atom  is  bigger.  This  may  cause  a  larger 
radiation  absorption  when  atoms  are  mostly  in  excited  states  as  in  the  case  of 


a  hot  plasma. 


oil 


Figure  2.  Elastic  cross  sections  for  hydrogen  a' 
| 200>,  |210>,  |211>,  |21-1>. 
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3.  L  N(a)  F  (q)  FOR  HYDROGEN  ATOM 
a^a  a' a 

As  mentioned  above,  F  ,  (q)  of  an  atom  or  the  scattering  cross 

section  of  electron-atom  collision  plays  an  important  role  in  the  energy 

intensity  emitted  due  to  electrons  moving  in  the  field  of  atoms.  Since 

we  shall  compare  the  radiation  absorption  coefficient,  which  is  related 

to  the  emitted  energy  intensity  (see  later),  to  the  measured  values 

in  a  hydrogen  plasma,  we  need  to  calculate  £  N(a)  F  ,  (q)  for  the 

hydrogen  atom.  Let  a  wave  function  of  a  hydrogen  atom  be  labeled  by  |nlm> 

where  n,i,  and  m  are,  respectively,  the  principal,  orbital  angular  momentum, 

and  magnetic  quantum  numbers.  Since  the  wave  function  of  a  hydrogen  atom 

2 

in  the  energy  level  En  has  n  degeneracies  for  a  spinless  orbital  electron, 
one  can  write 

Z  N(a)F  (q)  -  Z  ^  |<n '  I  'm*  |  l-elq*r| nim>  | 2  (IV-18) 

aa'  a  a  nira  2  4 
n'  /'ur  q 

where  K(n)  is  the  number  density  of  the  hydrogen  atoms  in  the  energy  E^ 

In  writing  eq.  (IV-18),  we  have  assumed  tyiat  the  states  of  the  hydrogen 
atoms  with  the  same  n  but  different  possible  values  of  l  and  m  are  equally  pop¬ 
ulated.  Expanding  the  absolute  square  in  eq.  (IV-18),  one  can  write 


Z  N(fi)F  (q)  =  —  Z^^  CG(n,rVl)  +  £  G(n,n/q)}  (IV-19) 

aa'  a  a  q  n  n2  n'yto 


where 

G(n,n,q) 


n-2Re  Z  <n/m|e  ^■~|nlm>  +  Z  |<nl'm'|e  ^*~|n/m>| 
im  Im 

I'm' 


(IV- 20) 
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G(n,n)q)  -  E  |<n'  fm'  |  ei‘2‘*-|nlm>  | 2.  (IV-21) 

tm 

r  m' 

Although  the  elastic  cross  section  for  hydrogen  atom  in  the  ground 

state  and  some  of  the  inelastic  cross  sections  when  the  atom  undergoes 

(17  22-27) 

certain  transitions  have  been  calculated,  *  there  are  no  explicit 

expressions  for  G(n,n,q)  and  G(n,n)q)  available  in  the  literature  except 

for  G( 1,1, q) ,  G(l,2,q),  G(l,3,q),  and  G(l,4,q).  In  order  to  get  G(n,n,q) 

and  G(n,n)q)  explicitly,  for  other  values  of  n  and  n'  we  shall  use  the 

(2S ) 

method  introduced  by  McCoyd,  Milford,  and  Wahl  which  we  present  below 
for  completeness. 

By  introducing  the  normalized  hydrogen  atom  wave  function  |n/m>  « 

N  nj(r)  Yjm(n)  and  expanding 

1 

ei9“-  -  E  [4n(2p+l)]2iPj_(qr)Y  (n), 

p-o  P  P° 

one  can  write 


<n'  / '  m'  I  e^  —  |nim>  = 


/+/' 


#m  #-m-  R „n#n.,.(<l)  (IV-22) 


where 


Y  In,  #•«.  “  iP[M2p+l)]2  !Y(Q)Y  (Q)Y*  (a)to 
p,  zm, /  m  po  zm  z  m  — 

■  (IV- 25) 


,1J2‘J5S 


with  (  being  the  Wigner  5J  symbol,  and 

m.  mm  1 
1  2 


(IV- 24) 
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p,ni,n'f"  ni,n'  l  'o  Jpw  n+f  n'+i'v 
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In  eq.  (IV-24),  j  (>x)  are  spherical  Bessel  functions.  L  are  the  associated 
p  n 


La  guerre  polynomials  and 
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ni  n' I ' 

o  o 
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(IV- 25) 


P  -  f+l'  +  2,  c'  =  2n/(n+n’),  c  ■  2n'/(n+n')  y  -  qaQnn '/(n+n' ) 
From  eqs.  (IV-22)  and  (IV-23),  one  can  perform  the  summations  over  m 


and  m'  in  eq.  (IV-21)  by  using  the  orthogonality  property  of  the  3j  symbol, 


(28) 


+£1  +^2  /W5\  /^1J2J3 ^ 

V"J1  ra2“’J2  \mim2y  V^lV^ 


2J3+1  Vi  VS' 


then  eq.  (IV-21)  becomes 


G(n,njq)  =  Z  |<n  '  I  'm'  |  e^  —  |  n/m>  | ' 
Im 
I'm ' 


(IV- 26) 

Pll\2 


Z  lRr,  ni  n  «  # '  (^)  I  (2p+l)(2/+l)  (2/  '+l)r  *) 

if  p  P>n*>n  *  OCO 

(29) 


When  the  Laguerre  functions  are  expressed  as  polynomials 


.21+1 


n-f-1 


.  .  ,  x  ^  .  vS+1  {(n+i^I) 

ntl  (0)‘io  (-1J  Tn:7thsT:-(9 


m2 


):(2i+i+s):s: 


ic  •  P 


(IV- 27) 
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the  radial  integral  in  eq.  (lV-2k)  reduces  to  sums  of  integrals  which  in¬ 
tegrate  directly  in  terms  of  hypergeometric  function^) 


_Y  P(^)P( P+P’+l) 

jx  e  j  (?x)dx  = 


,Pfl 


?E _  r/Trt'  +  l  P^l-P'  gP+3 

i(p+P’+l)  \  2  '  2  *  2  * 


(Kp*|)  d*,2)2 


The  hypergeometric  function  F(a,b,c,z)  is  the  analytic  solution  of  the 
hypergeometric  differential  equation ^ 

z( z-l)F'‘  +  [ (a+b+l)z-c]  F'+abF  =  o. 

About  the  singularity  z=o,  it  takes  the  form 


.  ab  a(a+l )b(b+l )  2  . 

r(a,b,o,7)  =  I  *  -  » *  2:c-('+i)  *  +  ••• 


-zi') 

l+72/- 


(IV- 2?) 


where  |z|<l-  Then  F(a,b,c«z)  will  be  a  polynomial  when  a  or  b  is  a  negative 

2/+1  ,  , 

integer.  For  the  terms  having  the  odd  power  of  x  in  the  product  of  L  ,  (ex; 

n+i 

and  (c'x),  one  can  prove  that  b  =•  is  a  negative  integer  and 

b<-l.  Therefore  the  radial  integral  of  eq.  (IV-28)  will  be  a  polynomial  for 
these  terms.  It  a  or  b  is  not  a  negative  integer,  F(a,b,c,z)  should  be 
expressed  about  the  singularity  at  z=l  instead  of  z=o 


F(u  b  c  y)  =  P( C ^  F(a  b  l>a+b-c  1-z) 
f  b>°>c>z'  f1(c-a|l(c-b)  na'D,i  8  D  c>x  z> 

(IV-50) 

+  Tn(iI)r^T'C)  ^1‘z^C*a"b  F(c-a,c-b,l+c-a-b,l-z) 
then  r’(a,b,c,z)  is  a  polynomial  in  (l-z)  because  c-a  is  a  negative  integer 

2J+1 

or  zero  for  the  terms  having  the  even  power  of  x  in  the  product  of  (cx) 
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2 1 '+]_ 

and  L  l+.,  (c'x).  The  first  term  in  eq.  (IV-30)  vanishes  because  p(c-a  )-*±» 
when  c-a  approaches  a  negative  integer  or  zero. 

With  the  tabulated  values  of  the  "3j"  symbols  and  the  above  described 
calculation  of  the  radial  integral,  one  can,  in  principle,  find  G(n,n,q) 
and  G(n,n^q)  for  any  values  of  n  and  n'.  It  is  obvious  that  the  calculations 
for  large  values  of  n  and  n'  are  tedious.  We  shall  calculate  G(n,n,q)  and 
G(n,n^q)  only  for  n,n'=l,2,3  below. 

From  eqs.  (IV-20),  (IV-22),  (IV-23),  and  (IV-26)  and  from  the  tabulated 

(32) 

values  of  the  "3j"  symbols,  one  can  obtain  (see  Appendix  D) 


G(2, 2, q)  = 


G(2,3,q)  = 


G(3,3,q)  - 


0(1,1, q)  =  U-Ro,io,io,(<i»2 

G(l,2,q)  =  R0,10,20^q  ^  +  5Rl,10,21^q^ 

G(i,3,q)  -  Ro,io,3o(q)  +  5Ri,io,3i(q)  +  5R§,  10,. 32^) 
**-2R0,20,20^q)  +  R0,20,20,(q)-6p0,21,21,(q' 

5R0,21,21^q^  +  6R2,21,21^q^  +  6R1, 20,21^ 

R0,20,30^q^  +  5R1, 20,31^  +  5R2, 20,32^  +  5R1, 21,30^ 

5R0, 21,31^  +  ^2,21, 31^  +  ^1,21,32^)  +  9R3,21,32(q) 

9-2R0,30,30^a)  +  R0,30,30^q^  _6R0,31,31^q^  +  5R0,3l,3l(q) 
6R2,3l,31(q)  "1OR0,32,32(q)  +  5R0,32,32(q)  +  ^  R2,32,32(q) 


(IV-31a 

(IV-31b) 

(iv-31c) 

(iv-3id) 

(lV-31e) 
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+  ^  R?,32,32^q^  +  ^1,30,31^  +  L0R2,30,!2^q)  +  12Rl,31,32^q) 


+  l8R3 ,31,32^* 


(IV-31f) 


where  H  (q)  are  determined  by  eqs.  (lV-2t)  thru  (IV-30).  From 

p,nf ,n  / 

straightforward  but  tedious  manipulations,  one  obtains  R  .  ,..(q)  for 

p,n/ ,n  x 


n  and  n'*l,2,3  with  the  following  results: 


Ro,io,io(q) 


R0,10,20(q) 


q  a 

(i  +  HT2) 
?  A* 


(9***q2a^)5 


Rl, 10,21^ 


R0, 10,30^ 


Rl, 10,31^ 


R2,10,32^q^ 


12  1 

2  2  qao 

2  5 

1 

32a2(^3t2) 


(^l2)4 

1  1  2 

,2  g2  S(^.L) 

'  2  1 

(^+|2) 

.2  ,2  .2 


51/2  <*|V 


R0,20,20^) 


nHx  2  2^ 
2q  a  -3q  a  +1 
o  o 

(l+q2a2)4 


Rl,20,2l(q) 


1  /  2  2  . . 
-  qa  (q  a  -1) 
2  o  o 

(i+q2a2)4 


(IV-32) 
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R0,21,21^ 


2  2 

q  a  -1 
o 

,,  2  2,k 

(i+q  80) 


R2,21,21^ 


„  2  2 
2q  a 

o 

,,  2  2,h 

(1+q  aQ) 


R0,20,30^q)  “ 


13  3 

2  232  2  75yU-  19i4y2+lL5 

56  7  (i+r2)5 


E1,20,31(<1)  6  (1+r2)5 


^2,20,32^  " 


aSif  72(L-:2l 
2  (in2)5 


Rl,2l,30^ 


12  r 

2  2  32  y(^7^-78 ?2+5) 
56  (1+72)5 


„10_2  2,,  2  v 

R  (q)  =  12,  7  122  =22 

ro,21,31U)  56  (1+r2)5 


11, 2  2f.  2  » 

R  (q)  =  2  2-  7  -I/- 

R2,21>51  q  ^  (l+72)5 


r  (q) ,  ££  7(?7.!-2i 

Rl>21,32  q  12  (1+  2}5 

5  2 

13, 2  5 

/  x  2  3  7 

R3,21,32Uj  "  5L3/2  (i+72}5 


(iv-32) 
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4  1+ 

/  \  2  5 

^Ji  32^^  s  "T”  o  £ 
4,52,52  3  (l+62)6 


3  6 

where  5  =  ^qaQ  f ?  *  5^ao*  '^le  substitution  e9s*  (IV-32)  into  eqs.  (lV-31) 
gives,  by  a  straightforward  manipulation, 


G(l,l,q)  «  {1-  ;  )2 


(IV-33) 


G(l,2,q) 


2i5q2a2 

(9+Uq2a2)5 


(IV-3I4) 


(b+xy 


(iv-35) 


G(2,2,q)  =  2x{—  + 


G(2, 


2  ( 

10 

5 

(l+x)2  ' 

a**)' 

(l+x)' 

375 

1160 

- -  +  - 

1728 

-G)  ,  x-q"aQ  (IV-36) 


5^  (l+x)'"’  (l+x)°  (l+x) 


G(3,3,q)  =  xC 


_9_  +  __2 _  +  _J _ +  83  _  236 


1+X  (l+x)2  (l+x)5  ( l+x ) ^  (l+x)5 


260  140  140  ,  922 

+  - r  +  - -  -  - 0}  ,  x  -  fq  ar 


(l+x)6  (l+x)T  (l+x)8'  ’  "  ^  ~° 

Equations  (IV-34)  and  (IV-35)  are  identical  with  the  results  obtained  by 

(22)  (17) 

R.  McCarroll  using  a  method  introduced  by  Mott  and  Massey.' 


(lV-38) 
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4.  RADIATION  ABSORPTION  IN  PARTIALLY  IONIZED  HYDROGEN  PLASMA 

With  expressions  of  G(n,n}q)  obtained  in  the  Last  section,  the  energy 

intensity  emitted  for  the  atom  undergoing  the  transition  from  the  level 

1  1  Eu 

En  to  Eni  and  for  Maxwellian  electron  distribution  2jt(jtQ)  NeEu  e  can  be 
obtained  from  eq.  (IV-7) 


Snn'(tt») 


Ne  ;V  T  fa.  /W  ^  2tgl  G(„,„;q) 


3N/2(n©m)5/2  6  0 


q  2 
n 


■  &(E  -E  ,+E  -E  ,-ftu) . 
n  n'  u  u' 


(IV- 39) 


Since  the  radiation  emission  coefficient  per  unit  length  for  an  isotropic 


medium  and  unpolarized  radiation  is  given  by 


E(to)  -  ^  /%  Z  E,(k), 


(lV-40) 


one  obtains  from  eq.  (IV-2) 


2  2 

®nn'(*°)  =  ^  Snn,^'ftD^ 


_  2L. 


•ft2n5ccc ? 


r  -  N  /d E  e‘  ®  jdE  ,  ;q"8X  53  »Sal  G(n,n;q) 

*  °  “  o  “  *d»  3 


■  6(En-En,+Eu-Eu,-ftJD) 


(IV-41) 


for  the  particular  emission  process  Eu+En-»Eu.+Eni-f‘fto.  Its  inverse  process, 
i.e.,  ^1+  En,+ftjj-*’Eu+En,  is  the  one  for  radiation  absorption,  and  the  absorp¬ 
tion  per  unit  length  can  be  obtained  by  interchanging  (u,u')  and  (n,n') 
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(te)  .  28.  ,  /aE  ,e"  6 

’  ^(nem)3/2*3  e°  *  ”  ”  5 


6<E„'-V  Eu’-Vto) 


(IV- 42) 


Then  the  absorption  ooefficient,  i.e.,  the  net  absorption  per  unit  length, 
is  (by  noting  that  G(n,n}q)*G(n',n,qy 


o'  .<*»)  -  A  .  (to)  -  E  „,<*»> 


JdE  /1e  ,  /q“‘^G(n,n'q) 
5^2  («9m)3/2“'3  e  0  u  0  U  q*i" 


Eu  Eu' -Eu 

N(n)  "  6  _n2  N(n')  ©  _n  B(E  .E  ,+E  -E  ,-ftu).  (IV-45) 

■  2  e  (  ,2  N(n)  6  J  V  "  n'  u  u 

n  n 

thermal  equilibrium,  the  number  densities  of  hydrogen  atoms  are  populated 


En' "En 

N(n')  _  Vf  e"  6 
N(n)  n2 


(IV- 44) 


a  ,(*)) 

nn 


7U) 

.  N  C  ( J-l)  /li  e"  9  /dEu.  C"  ^  ^  G(n’n;q) 
e  °  °  u  0  u  <W  q  n 


a  5(En-En.+Eu-Eu.-na>) 


(IV-45) 


where 


32  cV 


(n©m)  '  to' 


3/2  5  ‘ 


(IV-16) 
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Substituting  eqs.  (IV-33)  through  (IV-38)  into  eq.  (IV-15)>  one  finally 

1 

obtains  ann<(ha>)  as 


ann(fto)  =  Co(e°-l^hu  g^^Chu,©)  N(n)Ng 


(lV-i.7) 


for  n  =  n',  and  as 


a  ,(hx>)  =>  C  (e®-l)hi)  g  ,(hjo,0)  N(n)N 
nn'  o  Bnn  ’  e 


(IV- 18) 


for  nj=n'.  Here  the  upper  sign  of  gnnt  is  to  be  taken  when  n<n',  and  the  lower 
sign  when  n>n'.  The  definitions  of  gnn,  g^n'  and  gnn'  are  g^ven  below  with 
the  convention  that  g„n »  ~  gp'n  N(n')/N(n). 


2  ?  hu 

1  00  l+d+  6d_+9d_+l  6d^>9d++l  "  T  x 

g  (tin,9)  =  ~  /dx(  In  — “  +  - —  -  - ~ )e 

11  2  1  ld-  6(l+d.)5  6(Ud+)' 


(IV- 17a) 


1  „  l+b+  6b1+2lb5+16b2+Ub  +13 

g22(to,S)  -  ;  /dx{2In  —  *  - - 


13  2  hcu 

6b++2lb++16b++llb++13  -  —  x 

- 5 - )e 

3(i+b+r 


(IV- 17b) 


g53(*B'®)  =  Tq  {**  (9ini7r  + 


l+S+  5lS^+333S^+  1021S  i'+ 1150s  1056s  ?+6l5S_+ 177 


6(l+S.)7 


6  3  13  2  to 

51S++333S++1021S++1150S++1056S++615S++177  "  ©  x 

- - - }e 


gL2(to,0) 


6(l+S+) 


10  ±  ±  1  1 

?  V12  /idx{^ - - - i - 7 - i}e 

1  (9+l€'2b_)  (9+le'2b+) 


(IV-I7C) 


~~ x 


(!V-I8a) 


KMPWMnanraew 
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he\  S-+1*  3ef,s++^  *\  ‘  q~ 

+  81  ±  ±  r  13  _  .  __1 1 -  >e 

- t  l*\T±r^5  <> c 


13 


8  15  51  \k3s-+u>  (el/++4)- 


4(fta),0)  "  f23e23 


8  u  +  +  f375(ep,6.)-l78(e23^ 

^  M - ,  ±  ^6 

(1+e23P-} 


o 

5 


£”  flCD 

23 


375(6^+)-178(t^>599  ^  ^  — 


±  v6 

(1+e23P+) 


2ma0ficu 

5±  =  2~” 
ft 


d+sib±»  s±’?b' 


E  ,-E 

1  r 


P  +  55  ^  V  €nn'  =  •' J^T  ±l!  ’  fnn'  =  1 


f  ,  *  e 
nn' 


fro) 

6 


for  'flu)  <  E^ ,  -E^ 


E  ,-E 
n’  n 


for  fico>En,-En 


For*a»En,-En,  g^.fta,8)  “lu  be  evBlu8ted  by 

E  — E 

.  (f)9  6  £  (S-Ec^-c’^  VC)) 


E,-E, 


e'l5<'“"’9)  '  ,15 


=  £-  e  9  ^(12-6c2+Uc2  54-c2+c2eC2(6+c2)E1(c2) 


(IV- 18b) 


(IV- 18c) 


(iv-l7d) 


(IV- 18(1) 


(lV-U8f ) 


(IV-18g) 
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VE2 

AV-  9  J>(599  .Si 


-11 


•ft  03 


1283  2 

c_  +  — — —  c. 


5  3  30  3 


30  c3 


jQ6  I4  U8  5  A  b  c,,125  .  116  U8  2Nri  ,  \ % 

-  -g-  -  —  c'  +  c5  e  3(-^  +  __  y  .  c5)El(c5)) 


(lV-l(8h) 


where 


9ft 


2* 


23ft 


'1  2  7  P  2  f  ■*,  2 

32ma  0  9ma  0  ^  288ma  6 

o  o  o 


,  E^x)  =»  e"^  (iv-li8k) 


The  above  expressions  of  gj^t  and  g^n»  are  obtained  from  eq.  (rv-l+5)  con¬ 
sidering  the  cases  n<h’  and  n>n’  separately.  The  property  that  G(n,n,'q)=G(nJn,q) 
and  the  population  law  for  thermal  equilibrium,  i.e.,  eq.  (IV-l*4),  have  been 

used  in  obtaining  eq.  (IV-1j8).  The  values  of  E^(x)  are  available  in  tab- 

(3^  ± 

ulated  form.  .  g  .(ftcu.O)  and  g“  .(ftoo,©)  are  dimensionless  numbers  which 
nn  7  nn  7 

we  have  calculated  numerically.  Figures  3  through  11  show  the  results  of  the 
computations  for  g^n(ftci3,©)  and  ,  (ftoo,©)  as  functions  offta)  and  T.  g^,  is 
about  10  times  less  at  high  temperature  and  much  less  at  low  temperature  than 


5.  VALUES  OF  g  $103,0)  FOR  n>3. 
nn 

As  seen  in  sections  IV-3  and  IV-l»,  the  value  of  g  (fto3,0)  is  obtained 
through  the  calculation  of  G(n,n,q)  which  is  tedious  for  n>3.  Here  we 
shall  use  interpolation,  instead  of  the  direction  calculation  through  the 
method  in  section  IV-3,  to  get  the  values  of  for  n>3  by  knowing 

that  the  absorption  coefficient  due  to  inverse  bremsstrahlung  of  an  electron 
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in  the  field  of  a  neutral  atom  becomes  the  absorption  coefficient  due  to 
inverse  bremsstrahlung  cf  an  electron  in  the  field  of  an  ion  when  n  goes  to 
infinity,  i.e., 


lim 

n-*oo 


a  (*to,e)  bi,  . 
nn  ,  a  (<fau,Q) 


N(n)N 


NN 
I  e 

BI/ 


( IV- 1+9 ) 


where  N^.  is  the  number  density  of  ions  and  a  (ha),©)  is  the  absorption 
coefficient  due  to  inverse  bremsstrahlung  of  electrons  in  the  field  of  ions, 
has  been  given^  ^  for  hydrogen  atom  as 


BI,.  8  c  ,  ,3 

a  (fuu,e)  = - —  7-  (aaj 


E5/2 


NtN 


3(3*) 


1/2  A  KoJ  0l/2v3  Bff  "Ie 


(IV- 50) 


where  v  is  the  frequency  of  the  radiation,  E  is  the  ionization  potential  of 

a  hydrogen  atom  and  gff  is  the  free-free  Gaunt  factor v  *  '  depending  on 

temperature  and  the  absorbed  radiation  frequency.  With  the  value  of  gf^, 

one  can  find  the  upper  bound  of  g  when  n-*»  through  eqs.  (IV-U7),  (lV-49), 

nn 

and  (IV-50).  Since  we  shall  compare  our  calculated  absorption  coefficients 

with  the  measured  results  which  were  achieved  for  the  ruby  laser  frequency 

under  six  different  temperatures,  the  upper  bounds  are  obtained  for  these 

conditions  and  plotted  as  in  Figure  12  together  with  g  ,  g  ,  and  c 

n2  11  22  JJ 

obtained  from  Figures  3  through  5.  Then  the  values  of  g  for  n>3  can  be 

nn 

obtained  by  interpolation  on  the  smooth  curve  which  is  connected  through  the 


values  of  g^,  g^,  g^  ®nd  the  upperbound. 


Since  g  ,  is  about  ten  times  less  than  g^  when  n=l  and  is  getting 

much  less  than  g  when  n  and  n'  both  increase,  the  contributions  to  the 
nn 


g22  (Tiw.T) 


I04 


Figure 


(Tiw.T) 


-<h«*,T) 


I04 


Figure  7a.  g"2(£o,T). 


j 

I 


1 


(Ttw,T) 


-  (T»«,T) 


B  for  T*  2.2x10*  °K 
C  for  T*  2.0x10*  °K 
D  for  T *  1.5x10*  #K 


gnn(Hku,T)  vs.  l/n^  for  ruby  laser  frequency. 


81* 

absorption  due  to  g  for  n^n'  can  be  neglected  as  compared  to  the  contribu¬ 
tions  due  to  g  .  With  this  neglect,  the  absorption  coefficient  due  to  the 
nn  ’ 

sum  of  the  contributions  in  different  states  is 

a(fiD.e)  =  1  a  (fn>,9) 

n  nn 

ha) 

=  c  (e^  -l)fxa  N  l  g  ('faa,©)N(n)  (IV-51) 

o  e  ’  nn 

With  the  values  of  g  yui,©)  determined  above,  one  can  obtain  from  eq. 

nn  ’ 

(lV-5l)  the  radiation  absorption  coefficient  due  to  inverse  bremsstrahlung 

of  e  .ectrons  in  the  field  of  hydrogen  atoms  if  the  electron  density  N  , 

e 

the  neutral  densities  in  the  different  energy  levels  N(n)  and  the  temperature 
of  the  hydrogen  plasma  are  given. 

6.  EXPERIMENTAL  CONDITIONS  AND  THE  MEASURED  ABSORPTION  RESULTS 

Ir  .^66  Litvak  and  Edwards  ^ '  measured  the  absorption  coefficient  of 
the  ruby  laser  frequency  (X«691*3^)  for  different  initial  gas  pressures  in 
a  hydrogen  plasma  produced  by  a  giant  pulsed  laser  beam.  The  output  of  the 
laser  beam  for  producing  the  plasma  is  —  -  5  MW  peak  power  with  a  pulse 
width  of  about  18-3*  nseo.  A  25-cm  focal  length  lens  was  used  to  focus 
the  laser  output  near  the  center  of  a  brass  cubic  cell  which  contains  hy¬ 
drogen  gas.  The  initial  gas  pressures  before  the  -ight  went  in  were  It .7 , 

35 t  55 t  115  >  215,  and  1015  psi.  Although  the  initial  gas  temperature  was 
not  mentioned  in  Litvak  and  Edwards  work,  we  assume  that  it  was  room  tem¬ 


perature 
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The  measurements  of  the  absorption  coefficients  for  the  different 
initial  gas  pressures  were  performed  at  the  peak  luminosity  which  occurred 
near  the  end  of  the  laser  pulse.  The  electron  density  and  the  plasma  tem¬ 
perature  for  each  initial  pressure  corresponding  to  the  absorption  measure¬ 
ments  were  also  measured.  Table  I  shows  their  measured  results  of  plasma 

temperature  T,  electron  density  N  and  absorption  (aL)  .  for  different 

e  obs 

initial  gas  pressure  p^.  L  is  the  plasma  absorption  thickness  which  Litvak 
and  Edwards  assumed  to  be  varied  from  about  1cm  for  11.7  psi  to  1mm  for 
1015  psi. 

In  order  to  explain  the  measured  absorption  results,  Litvak  and  Edwards 
also  calculated  the  absorption  coefficient  0^  from  the  expression 


a  = 
v 


8 


3(3") 


_ i.  (CB  ,3 ill  N2(g 

1/2  fi  V  o  eVBff  ©  n  in 


E 

,n2© 


)(!-« 


fit!) 

0 


(IV-52) 


which  accounts  for  photoionization  and  inverse  bremsstrahlung  of  electrons 

( 54) 

in  the  field  of  ions.  Eq.  (IV-52)  is  obtained  from  the  absorption 
coefficient  of  a  hydrogen  atom  in  the  energy  level  E^  due  to  photoionization 


and  its  inverse 
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through  the  use  of  Saha  equation 


(IV- 53) 


N(n) 


f2nt£,3/2 
k  m©  ' 


2 

n 


exp 


\-QJ 


2 

e 


(lV-5»0 


06 


(  3^  35) 

where  g .  is  the  Gaunt  factor  for  free-bound  transitions/  *  rhe  cal- 
fn 

culated  a  shown  in  Table  I  is  two  orders  of  magnitude  less  than  the  mea- 
v 

sured  result  under  the  assumption  that  the  plasma  thickness  varies  from 
about  1cm  for  lb .7  psi  to  1mm  for  1015  psi.  This  large  descrepancy  in¬ 
dicates  that  the  measured  absorption  can  not  be  explained  by  the  photo¬ 
ionization  and  the  inverse  bremsstrahlung  of  electrons  in  the  field  of 
ions  through  the  use  of  Saha  equation. 

In  the  following  sections  we  propose  to  explain  the  absorption  which 
is  measured  in  this  experiment  by  considering  the  photoionization  process 
and  the  inverse  bremsstrahlung  of  electrons  in  the  field  of  neutral  atoms. 
In  this  explanation  we  shall  not  use  the  Saha  equation  to  predict  the  num¬ 
ber  of  neutral  atoms  in  the  plasma,  but  rather  we  shall  determine  it  from 
an  investigation  of  the  explosion  caused  by  the  laser  pulse. 


7.  DESCRIPTION  OF  AN  INTENSE  POINT  EXPLOSION 

It  was  shown  by  Litvak  and  Edwards  from  the  consideration  of  the  mea¬ 
sured  pressure  and  energy  variations  with  time  that  the  giant  laser  pulse 
produces  an  intense  p  int  explosion  with  a  spherical  shock  wave.  The  prob¬ 
lem  for  an  intense  point  explosion  has  been  investigated  by  Sedov, 
o 

After  the  energy  E  is  absorbed  into  the  gas  with  initial  pressure  p^  and 
mass  density  for  initiating  the  explosion,  a  shock  wave  forms  and  ex¬ 
pands  in  the  course  of  time.  Sedov  defined  p  ,  p  T  and  r  as  the  total 
mass  density,  the  pressure,  the  temperature  and  the  radius  of  a  point  be¬ 
hind  the  shock  wave  at  the  time  t  after  the  explosion.  Furthermore, 
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o  (Is  2 
r  (— ) 

Pi 


(IV- 55) 


and 


0  /E°\3 

r  =  (— ) 

Pi 


(iv-56) 


were  defined  and  Litvak  and  Edwards  called  them  as  the  characteristic  time 
and  shock  radius  at  which  the  counter  pressure  of  the  undisturbed  gas  nearly 
stops  the  expansion  of  a  spherical  explosion.  By  solving  the  equations  of 
motion,  continuity  and  energy,  one  can  determine  the  density,  pressure  and 

temperature  distributions  (—  and  —■)  as  functions  of  — ,  as  well  as 

?2  T2  r2  „ 

P2  Pg 

the  pressure  and  the  density  behind  the  shock  front  ( —  and  — - )  as  func- 

P1  pl 

r2 

tions  of  /  *  — .  The  symbols  p,p,T,  and  r  are  respectively  the  density, 
r° 

the  pressure,  the  temperature  and  the  radius  of  a  point  between  the  shock 
front  and  the  explosion  center  which  depend  on  time  implicitly  through 
r^(t).  Reference  (36)  contains  the  graphical  representations  of  the  above 
distributions  for  the  adiabatic  index  7=1. k.  For  our  later  use  and  for  a 
quantitative  understanding  of  their  variations,  we  reproduce  and 

p2  P2 

T 

—  in  Figures  13-15*  From  Fig.  13,  one  can  see  that  during  the  early  times 
T2 

after  explosion,  most  particles  are  concentrated  behind  the  shock  wave  and 
a  negligible  amount  of  particles  occupies  the  central  region  of  the  explo¬ 
sion. 

In  addition,  Sedov  also  obtained  the  following  equations 
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Figure  13-  Density  distribution  in 
a  point  explosion  behind  spherical 
shock  wave. 


Figure  l4.  Pressure  distribution  in 
a  point  explosion  behind  spherical 
shock  wave . 


Figure  15-  Temperature  distribution 
in  u  point  explosion  behind  spherical 

shock  wave. 
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for  spherical  shock  wave  where  q  is  the  square  of  the  ratio  of  the  sound 
speed  in  the  undisturbed  gas  to  the  shock  wave  velocity,  a  is  a  quantity 
which  depends  on  7  ana  the  shock  wave  geometry.  For  the  plasma  in  Litvak 
and  Edwards  experiment  7= 5/5,  and  for  a  spherical  shock  wave  one  can  find 
a  from  reference  (56)  as  0=2.  Then  one  obtains  from  the  above  mentioned 
equations 
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where  T»t/t°. 


By  taking  E°  to  be  the  energy  absorbed  from  the  laser,  Litvak  and 


Edwards  obtained  the  characteristic  shock  radius  r°  from  eq.  (IV-56)  which 
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is  given  in  Table  I.  The  initial  mass  density  0i  and  gas  pressure  pp  of 

the  plasma  are  known.  We  calculated  the  characteristic  time  t°  shown  in 

Table  I  from  eq.  (IV-55)  which  checks  with  the  Litvak  and  Edwards  result. 

In  addition,  Litvak  and  Edwards  determined  spectroscopically  the  peak 

luminous  volume.  The  absorption  was  measured  at  the  time  when  the  peak 

luminous  volume  occurred.  Although  the  luminous  volume  has  been  observed 

to  have  a  nonspherical  shape  due  to  the  rapid  axial  motion  occurring  during 

the  laser  absorption,  for  a  quantitative  discussion,  we  shall  assume  that  it 

has  spherical  shape  and  coincides  with  the  shock  volume  at  the  instant  when 

the  peak  luminosity  occurred  (we  shall  justify  this  assumption  presently). 

Then  the  radius  r^  )t  the  peak  luminous  volume  can  be  found  from  its  volume. 

The  values  of  r  f  r  each  initial  pressure  are  given  in  Table  I  together 

r2 

with  the  corresponding  values  of  l  =  —  .  Under  the  spherical  assumption 

r° 

of  the  peak  lumin  n;  lume,  r  turns  to  be  0.062cm,  much  less  than  1cm 
which  is  assumed  by  litvak  and  Edwards  as  the  absorption  length  in  the 
plasma.  We  shall  cal  late  the-  total  absorption  using  the  value  of  r^ 
as  obtained  a l  .ve  :r.  secti  n  -8. 

In  order  t  Justify  the  above  assumption,  we  now  calculate  the  time 
t  at  which  the  peak  luminosity  occurs  (i.e,  ,  the  time  at  which  the  absorp¬ 
tion  measurement  is  taken)  irom  eq.  (IV-57)  and  the  value  of  l  corresponding 
to  the  peak  luminosity  volume  in  Table  I.  The  results  are  also  given  in 
Table  I.  We  observe  that  t  varies  between  17-51  nsec.  As  we  mentioned 
before,  the  width  of  the  laser  pulse  varried  in  the  experiment  between 
18-36  nsec.  The  early  part  of  the  pulse  produces  the  plasma,  and  according 
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to  Litvak  and  Edwards,  the  peak  luminosity  occurs  near  the  end  of  the  pulse. 
Hence,  the  values  of  t  calculated  above  agree  reasonably  well  with  the 
experimental  conditions. 

Since  /  is  small  for  the  time  at  which  the  peak  luminosity  occurs, 
the  total  mass  density  o0  behind  the  shock  wave,  i.e.,  at  r0,  is  obtained 

<—  i— 

from  eq.  (IV-58)  to  be  about  four  times  the  initial  density  p^,  as  also 
pointed  out  by  Litvak  and  Edwards.  The  pressures  p^  behind  the  shock  wave 
at  the  time  the  peak  luminosity  occurs,  and  at  the  time  t=0.1psec  are 
obtained  from  eq.  (lV-59)  and  shown  in  Fig.  16.  For  comparison,  the  pressure 
measured  at  t=0.1psec  is  also  shown  in  the  same  figure.  The  measured  pres¬ 
sure  at  t=0.1(isec.  is  six  times  the  pressure  predicted  by  the  explosion 
theory  at  the  same  instant  for  the  initial  pressures  p^=  14 .7  and  35psi. 

At  higher  initial  pressures,  it  decreases  and  reaches  about  the  same  pres¬ 
sure  as  predicted  for  pj=1015  psi.  If  one  extrapolates  the  measured  pres¬ 
sures  up  to  the  time  at  which  the  peak  luminosity  occurs,  a  descrepancy  is 
observed  between  the  extrapolated  value  and  the  presssure  predicted  by  the 
point  explosion  theory  for  the  same  instant.  Here  also,  this  descrepancy 
is  large  for  low  initial  gas  pressures  and  small  ior  high  initial  pressures. 
This  may  indicate  that  the  point  explosion  assumption  is  better  justified 
at  high  initial  gas  pressures  than  at  low  initial  pressures.  Since  the 
focal  volume  of  the  laser  beam  is  elongated  in  the  direction  of  the  beam, 
and  since  the  energy  absorption  decreases  away  from  the  source  due  to  the 
attenuation  of  initial  laser  beam  producing  the  shock,  an  egg-shaped  shock 


front  is  perhaps  a  more  accurate  description  than  the  spherical  shock  front 
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as  we  assumed  above.  Furthermore,  the  initial  volume  of  explosion  which  can 

"5  3 

be  taken  as  the  focal  volume  (10  cm')  of  the  laser  beam  is  not  negligible  as 
compared  to  the  shock  volume  ( 8  x  10  ^cm  even  at  0.1  nsec.  This  may  also  be 
a  contributing  factor  to  the  discrepancy  in  the  measured  and  calculated  pres¬ 
sures  . 

The  point  explosion  assumption,  as  pointed  above,  predicts  the  lo¬ 
cation  of  the  shock  front  reasonably  well.  Since  we  need  the  result  of  the 
point  explosion  theory  only  to  estimate  the  neutral  density  distribution 
within  the  shock  volume  and  the  shock  dze  in  our  absorption  calculations, 
and  since  the  pressure  behind  the  shock  front  does  not  enter  our  absorption 
formulas,  the  above  discrepancy  in  the  estimation  of  pressures  is  not 
critical  for  our  purpose. 

8.  ABSORPTION  CALCULATION  BASED  ON  NEUTRAL  DISTRIBUTION  INSIDE  A  SHOCK 
WAVE 

In  this  section,  we  shall  calculate  the  absorption  due  to  the  photo¬ 
ionization  process  and  the  inverse  bremsstrahlung  of  electrons  in  the  field 
of  neutral  atoms  by  using  the  particle  distribution  inside  a  shock  wave, 
i.e.,  Fig.  13.  In  doing  this,  one  has  to  determine  first  the  population 
of  the  excited  levels  of  the  hydrogen  atoms  as  a  function  of  position  in 
the  plasma.  We  assume  local  thermal  equilibrium  among  the  neutral  atoms, 
i.e.,  eq.  (iV-UU)  so  that  the  relative  populations  of  the  excited  levels 
is  not  an  explicit  function  of  position  (it  may  depend  on  position  through 
temperature).  Then  the  local  absorption  coefficient  due  to  the  inverse 
bremsstrahlung  of  electrons  in  the  field  of  neutrals  can  be  written 
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aBN(r)  =  An(a o,6)  Ng(r)  K(l,r)  (IV-60) 

2  -E  1 0 

where  N(n)/N(l)  =  n  e  n  and  A  (oj,0)  is  a  function  of  temperature  and 

radiation  frequency.  An(o>,9)  increases  with  n,  approaching  finite  value 

A^(a>,©)  as  n-*tc.  Since  En-*E  as  n-*»,  the  above  summation  diverges  unless  it 

is  truncated  at  some  n=n*.  The  physical  reason  for  truncating  the  summation 

at  some  n=n*  can  be  explained  as  follows:  Because  of  the  interaction  of 

an  atom  with  the  nearby  particles  in  a  plasma,  the  ionization  potential 

will  be  lowered  when  the  atom  is  inside  a  plasma.  An  excited  atom  in  a 

level  above  n*  must  be  treated  as  an  ion  and  a  free  electron  even  though 

the  level  may  be  below  the  ionization  potential  of  an  unperturbed  atom. 

Several  theories  for  the  determination  of  n*  and  the  lowering  of  the 

ionization  potential  AE  have  been  proposed  in  literature.  Drawin  and 
(37 ) 

Felenbok  have  reviewed  these  theories  and  showed  that  they  all  yield 
similar  results,  i.e., 

2  Pn 

n*  *  (S+l)  “  (IV-61) 

ao 


AE  ^  (S+l)  _£. 

PD 


where 


is  the  Debye  radius 


(  38,  39) 


given  by 


PD  = 


(IV- 62) 


(IV- 63) 
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In  eqs.  (lV-6l)  and  (TV-62),  S  is  the  ionization  stage  of  the  particle 
under  consideration  and  S=0  for  neutral  atoms.  The  quantities  and 
in  eq.  (IV-63)  are  respectively  the  density  and  charge  of  neighboring  ions  and 
for  hydrogen  plasma,  z^=l  and  N^=Ng.  We  have  calculated  and  tabulated  in 
Table  I  the  values  of  n*  for  each  initial  gas  pressure  using  the  measured 
electron  density  and  temperature  in  the  experiment  by  Litvak  and  Edwards. 

After  having  determined  n*  in  eq.  (IV-60),  we  now  consider  the  posi¬ 
tion  dependence  of  the  number  density  in  eq.  (IV-60).  According  to  the 
intense  point  explosion  theory,  the  relative  mass  density  distribution 
follows  the  curves  in  Fig.  13  for  the  early  times  after  explosion.  One 
observes  that  the  particles  occupy  vitually  a  very  narrow  spherical  shell 
of  a  thickness  of  the  order  of  0.2  rg.  We  shall  refer  to  this  region  as 
the  shell  region  in  below.  The  central  region  contains  very  few  particles. 
However,  the  temperature  increases  very  rapidly  in  this  region  towards 
the  center  and  the  pressure  is  constant  with  a  value  of  approximately 
We  shall  refer  to  this  region  as  the  central  hot  core.  Let  the  relative 
number  density  of  the  atom  at  the  point  r  between  the  shock  wave  and  the 

explosion  center  be  cp(r )*N(r)/N,p  and  the  relative  electron  density  at  the 

n* 

point  r  be  9~(r)=N  (r)/N  where  N(r)  *  T  N(n,r)  and  NP  and  Ne  are  respec- 
e  e  n=l  eL 

tively  the  total  neutral  density  behind  the  shock  front  and  the  average 
electron  density  inside  the  shock  volume.  From  the  conservation  of  the 
particles 

r2  .  , 

Ijn  / drr  [N^cp(r)+Necpe(r)]  *  ~  r^N^ 


where  H  is  the  t  *  .  nur.  1  er  aensity  of  hydrogen  atoms  before  the  explosion 
(i,e.,  before  the  r.ciaence  )f  the  laser  beam)  and  r  is  the  radius  of  the 
shock  wave.  After  hanging  the  variable,  we  have 

'2  Ni  Ne  1  2 

/cixx  vp(x)  =  —  -  —  / dxx  q>  (x)  ( IV-61* ) 

3^  n2  0  e 

and 


where 


.1  p 

3  j  dxx  9  ( x )  =  1 
0  e 


(IV-65) 


X 


-5  3 

The  laser  beu  .d.;  i  >aused  within  a  volume  of  10  cm  which  is  smaller 

-3  3 

than  the  peak  lutr  m*  volume  of  1  cm  at  the  time  of  the  absorption 

measurement.  The  r  .  of  the  laser  beam  focal  volume  is  about  5  times 

less  than  the  rant  ' T’  the  peak  luminous  volume.  No  geometry  correction 
is  needed  in  defer'  z  the  optical  path,  hence  the  absorption  due  to 
the  inverse  hr-  it :•  g  f  electrons  in  the  field  of  neutral  atoms  can 
be  expressed,  r  ir  ;.cing  the  relative  number  densities  into  eq.  (lV-60), 
as 


V  /dxAn(ui,0)N  W  r  cp  (x)cp(x).  ( IV-66 ) 

---1-1  e  2  N  2  e 

Equation  (iV-co  ar.aot  be  easily  calculated  without  any  assumption 
about  q?e(x),  c(x  *  '  emperature  distribution.  The  problem  will  be 

complicated  when  t  e  t  mreratare  distribution  (see  Fig.  15)  in  the  shock 


97 


wave  as  a  function  of  position  is  considered.  The  plasma  temperature  and 
the  electron  density  are  measured  through  the  line  broadening  of  H  which 
depends  upon  the  electron  and  neutral  densities.  Since  most  of  the  hydrogen 
atoms  are  confined  in  the  shell  region  behind  the  shock  front,  the  electron 
density  and  the  plasma  temperature  measured  through  the  line  broadening  is 
more  likely  to  indicate  the  average  temperature  and  electron  density  in 
the  shell  region.  Furthermore,  the  degree  of  ionization  is  not  uniform 
in  the  shock  volume.  Due  to  the  very  high  temperatures,  the  hydrogen  gas 
can  be  expected  to  be  fully  ionized  in  the  hot  core.  However,  the  degree 
of  ionization  is  more  likely  to  decrease  towards  the  shock  front  because 
the  temperature  there  is  of  the  order  of  10  °K  and  thus  not  sufficient  for 
onization.  One  may  conclude  from  this  argument  that  the  electron  density 
wi.l  also  be  a  decreasing  function  of  the  radius.  However,  the  rapid  in¬ 
crease  in  the  particle  density  towards  the  shock  front  may  result  in  a 
uniform  electron  density  distribution  in  spite  of  the  decrease  in  the  de¬ 
gree  of  ionization.  Hence  we  shall  assume  that  the  temperature  and  the 
electron  density  are  uniform  in  the  shell  region  with  the  measured  values. 

The  contribution  to  the  absorption  in  the  hot  core  where  we  expect 
almost  full  ionization  is  negligibly  small  as  compared  to  the  measured 
absorption.  The  mechanisms  responsible  for  photon  absorption  in  this 
region  are  the  inverse  bremsstrahlung  of  electrons  in  the  field  of  ions  and  the 
photoionization.  In  fact,  since  the  gas  is  almost  fully  ionized,  the  relative 
contribution  of  the  photoionization  process  is  small  as  compared  to  the  in¬ 


verse  bremsstrahlung  of  electrons  in  the  field  of  ions.  Litvak  and  Edwards  calcu- 
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lated  the  absorption  in  this  region  using  the  measured  electron  density  and 
temperature,  and  assuming  the  Saha  equation  in  estimating  the  neutral  density. 
Even  if  one  assumes  a  plasma  thickness  of  1cm,  their  result  account ed  only  one 
percent  of  the  measured  absorption.  But  according  to  explosion  theory  the 
size  of  the  hot  core  is  of  the  order  of  0.12cm.  so  that  the  absorption  is 
less  than  0.1$  of  the  measured  value. 

We  thus  calculate  that  the  main  absorption  takes  place  in  the  shell 
region  where  the  gas  is  partially  ionized.  Under  these  assumptions,  eq. 

(IV- 66)  becomes 


(aL)™=  Ta  (»,e)N  H  ^  r  /dxtp(x), 
n=l  n  e  d  H  ^  -1 


(iv-67) 


where  we  include  the  hot  core  also  for  convenience  even  if  i*.s  contribution 
is  small.  The  neutral  density  distribution  cp(x)  is  shown  in  Fig  13  and 
has  an  appreciable  value  only  for  x  close  to  1.  It  is  sufficient  for  our 
purpose,  but  not  necessary,  to  approximate 


/"  cp(x)dx  /”cp(x)x  dx  =  — 7; 


N,  -N 


3N„ 


Then  eqs .  (IV-67)  reduces,  in  view  of  eq.  (lV-51),  to 


/  t xBN  2  x 

(ctL)  =  -  c0hjo 


/—  \  n* 

(XV-68) 


where 


=  fznV1^9 


n=l 


is  the  truncated  partition  function  of  hydrogen  atom  with  the  calculated 


J 
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TABLE  I 

MEASURED  AND  CALCULATED  ABSORPTIONS  FOR  RUBY 
LASER  FREQUENCY  IN  THE  HYDROGEN  PLASMA 


PL(psi) 

14.7 

35 

55 

115 

215 

1015 

T(°K) 

2.2x10? 

2.0x10? 

2.1*xl05 

1.5xl05 

l.l*xl05 

1 

9.0x10 

N  (cm,^) 
e 

1+.5x101R 

1019 

5xl018 

1019 

8.5xl018 

1019 

t°(nsec. ) 

b.9 

4.3 

4.3 

1*.0 

3.2 

1.7 

r°(cm. ) 

0.52 

0.1*6 

0.1*6 

0.1*3 

O.31* 

0.18 

r2(cm. ) 

0.062 

0.062 

0.062 

0.062 

0.062 

0.01*2 

l 

0.12 

0.13 

0.13 

o.ll* 

0.18 

0.23 

t(nsec. ) 

17 

18 

18 

21 

31 

27 

n* 

11* 

12 

11* 

11 

11 

10 

N  -N  (cmT3) 
1  e 

l*.9xl019 

1.2xl0?° 

20 

2.0x10 

l*.lxl02° 

7.7xl020 

3 .7x10 

z 

998 

598 

1062 

362 

337 

11:0 

a  (cmL) 

V 

0.012 

.0.011 

0.020 

0.15 

0.12 

0.6 

(aL)obs. 

0.72 

0.80 

1.9!* 

3-0 

2.61* 

1*.  12 

(ol)bn 

0.015 

0.087 

0.062 

0.1*7 

0.75 

*9 

(0L)PI 

0.01*7 

0.20 

0.18 

1.13 

2.26 

11*. 9 

(ol;BN+PI 

0.062 

0.29 

0.21* 

I.60 

3.01 

19.8 
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values  given  in  Table  I. 

In  a  similar  way,  one  can  obtain  the  absorption  due  to  the  photoioni¬ 
zation  and  its  inverse  from  eq.  (IV-53)  as 


fijo 

z.xPI  128  2  ,E^.  ‘  Tx  VNe  ^*2gfn  -En/0 

<aL)  '  ^7i  “o  (ta>  (1-S  >  r2  —  nil—  e  .  dv.69) 

Taking  r^  as  the  radius  of  the  peak  luminous  volume  which  is  assumed  to 

BN  PI 

be  spherical,  we  have  calculated  the  absorptions  (aL)  and  (OL)  from 


eqs.  (IV-68)  and  (IV-69).  The  results  are  given  in  Table  I  together  with 

..  .  ,  T  xBN+PI 

their  sum  (aL) 


9.  DISCUSSIONS 

BN+PI 

From  Table  I,  one  can  see  that  the  calculated  absorptions  (aL) 

due  to  the  photoionization  and  the  inverse  bremsst rah lung  of  electrons  in 

the  field  of  neutron  atoms  are  not  always  in  good  agreement  with  the  measured 

result  (aL)  ,  .  It  increases  from  the  value  of  ten  times  less  for  P ,  =  lU .T  psi 

obs  1 

to  the  value  of  five  times  larger  for  p^=1015  psi  than  the  measured  absorp¬ 
tion.  According  to  the  description  given  by  Litvak  and  Edwards,  the  shape 
of  the  luminous  volume  was  not  exactly  spherical.  Although  their  description 
was  not  explicit  enough  for  their  experiment,  they  referred  to  other  similar 
experiments  in  which  the  cigar-shaped  or  egg-shaped  luminous  regions  were 
observed  in  the  direction  of  the  laser  beam.  Furthermore,  the  discrepancy 
between  the  measured  pressure  and  that  predicted  by  the  intense  point  ex¬ 
plosion  theory  with  a  srsherical  shock  wave  is  larger  at  low  initial  gas 
pressure  than  at  high  pressure.  This  may  suggest  that  the  shape  of  the  shock 
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volume  at  low  initial  gas  pressure  deviates  more  from  spherical  at  high 
initial  pressure.  If  we  intepret  r^  as  the  major  radius  of  the  luminous 
volume,  then  we  may  predict  a  larger  absorption  than  we  calculated  by  using 
a  spherical  luminous  volume.  In  fact,  this  may  be  the  reason  why  Litvak 
and  Edwards  assumed  the  absorption  thickness  of  1cm  for  p^=1^.7  psi.  If 
this  assumed  absorption  thickness  is  correct,  r^  would  be  0.5cm,  instead 
of  0.062cm,  and  the  calculated  absorption  with  r^=0.5cm  will  be  almost  the 
same  as  the  measured  value  at  p^=ll+.7  psi,  provided  that  the  neutral  density 
distribution  along  the  major  axis  has  a  similar  distribution  to  that  in  the 
spherical  case  shown  in  Fig.  15.  At  any  rate,  our  interpretation  predicts 
the  absorption  better  than  a  factor  10,  in  fact  in  most  cases  even  better 
than  a  factor  of  6.  Furthermore,  our  calculation  explains  the  increase  of 
the  absorption  with  the  initial  pressure,  independently  of  the  possible 
dependence  on  the  pressure  of  the  apparent  plasma  size. 

Another  factor  which  may  be  contributing  to  the  above  discrepancy  is 
the  assumption  of  a  uniform  temperature  and  electron  density  distribution 
in  the  shell  region  behind  the  shock  front.  However,  the  error  due  to 
this  assumption  is  not  expected  to  be  significant,  because  we  have  found 
only  a  decrease  by  a  factor  of  2/5  assuming  a  linear  electron 
density  distribution  in  the  shell  region  and  taking  the  electron  density 
to  be  zero  at  -the  shock  front. 

The  calculated  absorption  due  to  the  photoionization  and  its  inverse 
is  about  three  times  the  absorption  due  to  the  inverse  bremsstrahlung  of 
electrons  in  the  field  of  neutral  atoms  for  n*  =  10.  It  depends  on  n*. 
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For  small  n*,  its  contribution  is  dominant  and  for  large  n*  it  is  negligible 
as  compared  to  the  inverse  bremsstrahlung  of  electrons  in  the  field  of  neutral 
atoms.  Furthermore,  photoionization  can  occur  only  for  the  levels  above  and 
equal  to  n  such  that  E-I^<  ftjo  and  the  absorption  due  to  this  process  varies 
as  l/n5 •  For  low  energy  photons,  such  as  the  carbon  dioxide  laser  (XalO^xlO*4 
cm,  ftu=0.117ev),  its  contribution  will  be  negligibly  small  as  compared  the 
absorption  due  to  the  inverse  bremsstrahlung  of  electrons  in  the  field  of 
neutral  atoms. 

The  absorption  per  electron  and  per  neutral  atom  in  any  level  n  due  to 
the  inverse  bremsstrahlung  of  electrons  in  the  neutral  atoms  is  smaller  than 
the  absorption  per  electron  and  per  ion  due  to  the  inverse  bremsstrahlung 
of  electrons  in  the  ions  as  indicated  in  Fig.  12.  They  are  of  the  same 
order  of  magnitude  for  almost  all  n  and  becoming  equal  when  n-*».  The  rel¬ 
ative  importance  of  these  two  mechanisms  depends  on  the  ratio  of  the  neutral 
to  ion  densities  N  /N. .  If  N  /N,  is  a  large  number,  the  inverse  bremss- 
trahlung  of  electrons  in  the  field  of  neutral  atoms  is  important,  otherwise 


it  is  small. 


CHAPTER  V 


CONCLUDING  REMARKS 


In  the  first  part  of  this  thesis  we  have  compared  the  photon  transport 
theory  in  dispersive  media  to  the  Maxwell's  wave  theory  by  considering  the 
index  of  refraction  and  the  photon  absorption  per  unit  time  as  well  as  per 
unit  length.  In  the  photon  transport  theory  the  effect  of  the  medium  is 
taken  into  account  by  assigning  a  different  frequency  to  photons  of  a  given 
wave  number  in  the  medium,  than  their  frequency  in  vacuum.  It  is  implied  in 
this  theory  that  the  wave  number  is  the  same  in  the  medium  and  in  vacuum. 

In  Maxwell's  theory,  the  effect  of  the  medium  is  characterized  by  a 
functional  which  relates  the  macroscopic  current  to  the  electric  field. 

When  linearized,  this  functional  is  completely  defined  by  the  conductivity 
tensor  in  the  transformed  (k,s)  domain.  The  conductivity  tensor  is  obtained 
quite  generally  by  using  Kubo '  s  linear  response  theory  in  terms  of  the  micro¬ 
scopic  currents.  Thus,  we  can  calculate  the  index  of  refraction  and  the 
damping  coefficient  both  in  time  and  space  in  the  framework  of  the  Maxwell's 
theory  first  in  terms  of  conductivity  and  then  in  terms  of  the  microscopic 
currents  with  Kubo's  theory.  In  other  words,  we  can  express  the  above 
observable  macroscopic  quantities,  in  terms  of  microscopic  quantities  through 
the  Maxwell's  equations  which  describe  the  electromagnetic  phenomena  in 
arbitrary  media  macros copically. 

It  is  at  this  stage  one  can  compare  the  photon  transport  theory  to  the 
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Maxwell's  wave  theory,  because  in  the  former  the  index  of  refraction  and 
the  damping  coefficient  both  in  time  and  space  are  expressed  in  terms  of  the 
microscopic  currents. 

Following  the  above  procedure,  we  have  found  that  both  theories  yield 
the  same  results  for  the  index  of  refraction  and  the  damping  coefficient  per 
unit  time  only  in  the  weakly  absorbing  media.  When  the  medium  is  strongly 
absorbing,  the  results  look  quite  different  although  we  have  not  estimated 
the  difference  numerically  in  specific  problems. 

As  to  the  damping  in  space,  the  expressions  obtained  from  the  two  theories 

for  the  index  of  refraction  and  the  damping  coefficient  per  unit  length  are 

similar  only  if  the  medium  is  both  weakly  absorbing  and  slightly  dispersive. 

More  explicitly,  the  damping  coefficient  per  unit  length  is  obtained  in  a  weakly 

absorbing  medium  (i.e.,  n,  «  n„)  as  — —  af(n_k,u»)  and  — —  a^(k,  — )respec- 

’  1  o  n  c  \  0  ’  n„c  \  ’  n_ 

o  o  o 

tively  in  the  Maxwell's  theory  and  the  transport  theory.  Clearly  if  n^l, 
the  results  are  identical. 

We  feel  that  a  better  correspondence  between  the  transport  and  wave 
approaches  in  regard  to  the  damping  in  space  can  be  established  if  the  photons 
are  dressed  such  that  the  frequency  is  required  to  be  he  same  but  the  wave 
number  is  allowed  to  be  different  in  the  medium.  More  research  in  this 
direction  seems  to  be  called  for. 

The  radiation  absorption  due  to  inverse  bremsstrahlung  of  electrons  in 
the  field  of  neutral  atoms  is  formulated  by  using  free  electron  wave  functions. 
In  this  approximation,  the  calculations  involving  atom-electron  cross  sections 
turn  out  to  be  identical  to  the  use  of  the  first  Born  approximation.  The 
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elastic  cross  section  for  hydrogen  atoms  in  the  ground  state  calculated  in 
Born  approximation  is  less  than  the  cross  section  calculated  from  the  partical 
wave  method.  This  discrepancy  is  large  at  specially  low  electron  energies. 
Since  the  elastic  cross  section  predicted  by  the  partial  wave  method  is  in 
good  agreement  with  the  experiment,  a  more  accurate  formulism  of  the  absorp¬ 
tion  problem  can  be  achieved  by  using  the  partial  wave  method  rather  than 
plane  waves.  However,  since  the  elastic  and  inelastic  collisions  of  electrons 
with  the  atoms  in  excited  states  are  also  involved  in  the  problem,  the  use 
of  the  partial  wave  method  would  make  the  problem  much  too  complicated. 
Undoubtedly,  the  use  of  free  electron  wave  functions  in  the  problem  will  intro¬ 
duce  some  error,  perhaps  predicting  smaller  values  for  the  absorption.  How¬ 
ever,  at  this  stage,  one  is  satisfied  with  an  order  of  magnitude  agreement 
between  the  measured  and  calculated  absorptions  due  to  the  uncertainties  of 
the  experimental  condition.  This  justifies  the  approach  we  have  taken  in 
this  work. 

In  the  formulation  of  the  radiation  absorption,  the  second  quantization 
is  used  to  express  the  potential  between  atom  and  electron  as  well  as  the 
interaction  between  particle  and  radiation  in  terms  of  the  particle  and 
radiation  creation  and  destruction  operators.  In  this  process  we  considered 
only  the  binary  collisions.  In  addition,  we  assumed  the  atoms  in  the  medium 
to  be  infinite  heavy.  These  assumptions  are  adequate  for  our  stated  purpose 
in  this  thesis. 

The  neutral  density  determined  through  the  Saha  equation  (which  holus 
when  all  the  particles  are  in  local  thermal  equilibrium)  and  using  the 
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measured  electron  density  and  temperature  in  Litvak  and  Edwards  experiment  is 

about  three  orders  of  magnitude  less  than  t'he  initial  particle  density  The 
absorption  calculated  by  Litvak  and  Edwards  in  considering  only  the  photo¬ 
ionization  and  the  inverse  bremsstrahlung  of  electrons  in  the  field  of  ions 
through  the  use  of  Saha  equation  is  negligible  as  compared  to  the  measured 
result . 

According  to  intense  point  explosion  theory,  most  particles  are  concen¬ 
trated  in  the  shell  region  behind  the  shock  front  and  a  negligible  amount 
of  particles  occupies  the  central  hot  core.  Local  thermal  equilibrium 
among  the  neutral  atoms,  instead  of  Saha  equation,  is  assumed  in  this  shell 
region.  The  absorptions  calculated  by  considering  the  photoionization  and 
the  inverse  bremsstrahlung  of  electrons  in  the  field  of  neutral  atoms  in 
this  shell  region  are  not  in  very  good  agreement  with  the  measured  results. 
However,  the  agreement  is  always  better  than  a  factor  of  10  and  in  fact 
better  than  a  factor  of  6  in  all  but  one  case.  The  discrepancies  between 
the  calculated  and  measured  results  may  be  attributed  mainly  to  the  use 
of  the  radius  of  the  peak  luminous  volume,  which  is  assumed  to  be  spherical, 
as  the  actual  shock  radius. 

In  this  thesis  we  have  calculated  the  absorption  coefficient  per  unit 
atom  and  electron  due  to  the  inverse  bremsstrahlung  of  electrons  in  the 
field  of  neutron  atoms  as  a  function  of  the  electron  temperature  and  radia¬ 
tion  frequency,  and  presented  the  results  graphically.  With  these  curves 
and  the  conventional  formula  for  photoionization,  one  can  now  estimate  the 


total  absorption  due  to  the  photoionization  and  the  inverse  bremsstrahlung 
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of  electrons  in  the  field  of  neutrals  if  the  neutral  density  in  the  plasma 
is  Known. 

The  interpretation  of  the  Litvak  and  Edwards  absorption  measurements  in 
this  thesis  by  considering  only  the  above  absorption  mechanisms  due  to  the 
neutral  atoms  and  using  the  point  explosion  theory  with  certain  plausible 
arguments  to  guess  the  electron  density  and  temperature  distributions  is  mainly 
suggestive.  In  the  absence  of  any  accurate  information  for  the  plasma  size  and 
of  an  explosion  theory  which  takes  into  account  the  finite  initial  volume  of 
the  explosion,  the  agreement  obtained  in  this  work  between  the  measured  and 
calculated  absorptions  is  considered  as  a  strong  evidence  for  the  importance 
of  neutral  atoms  in  the  interpretation  of  the  absorption  experiments  in  plasmas. 
In  fact,  the  absorption  due  to  the  neutrals  may  even  be  the  dominant  heating 
mechanism  which  causes  the  explosion. 

More  experimental  work  designed  primarily  for  the  verification  of  the 
importance  of  the  neutral  atoms  is  needed. 


APPENDIX  A 


DERIVATION  OF  EQUATION  (II-U2) 


In  Chapter  II,  we  already  obtained  the  following  expressions 
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The  substitution  of  eq.  ( 11-29)  into  the  above  expressions  gives  after  using 
eqs .  (II-23)  and  ( 11-30), 
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Equation  (A-4)  and  the  first  sum  in  the  bracket  of  eq.  (A-3 )  comes,  respectively 
from  letting  n  =  1^  in  eq.  (A-2)  and  n'=q+l  in  eq.  (A-l).  The  second 

the  bracket  of  eq.  (A-3)  comes  from  letting  V«t,-1  in  e**  (A_l)  *  Eva1uating 

the  matrix  elements  of  A(r)  and  A(r,t)  by  eqs .  (11-12)  and  (H-33)  and  assuming 

5 

that  the  radiation  field  does  not  change  appreciably  over  L  ,  then 
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Interchanging (r,r')  and  letting  t'=-t'  in  the  second  term  of  eq.  (A-6),  and 
summing  the  intermediate  states  I n *  >  and  |n">,  one  finally  obtains  eq.  (11-42) 
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APPENDIX  B 


SECOND  QUANTIZATION 

In  this  appendix,  we  shall  express  the  potentials  Vle,  Z1  and  the 

interactions  V.  ,  V  ,  and  V.  in  terms  of  the  particle  creation  and  destruction 
operators  by  second  quantization.  In  section  III-l,  V^e  represents  the  sum 
of  the  potentials  between  an  electron  and  an  atom,  i.e., 

Ze  =  rZf 

am 

pm 

where 
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is  the  potential  between  the  o-th  atom  and  the  m-th  free  electron  with 

and  RCT  being  the  positions  of  the  free  electron  and  the  nucleus.  p^  is 
the  position  of  the  j-th  atomic  electron  with  respect  to  R_.  In  second 
quantization,  can  be  written  as 
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where  A(u)  is  the  destruction  operator  of  destructing  a  free  electron  of 
momentum 'flu  and  A(Ka)  is  that  of  destructing  an  atom  of  the  external  momentum 
'fiK  and  the  internal  state  *a>.  A^(u')  and  A^(K'a')  are,  respectively,  the 

creation  operators  of  creating  a  free  electron  of  momentum 'flu '  and  an  atom 
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of  external  momentum  ^K'  and  the  internal  state  ^a'> 


The  wave  function  of  a  free  electron  is  given  by 
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(B-2) 


and  the  wave  function  of  describing  the  external  and  internal  states  of  an 
atom  is  given  by 
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where  we  have  assumed  that  the  center-of-mass  coordinates  ol'  the  atom  coin¬ 
cide  with  its  nucleus  coordinates.  Substituting  eqs.  (B-2)  and  (B-3)  into 
eq.  (B-l),  one  obtains 
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where  the  subscript  K  of  denotes  the  Kronecker  delta, 


In  the  same  way,  one  can  obtain 
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£a£'a'  L5  | K- K ’  |  2  J=1 

/b/'b' 


<b  ’  ( z-  I  e^- —  |b  >A+(K’a'  )A+(_/'b'  )Av'Ka)A(_ib)  (III-3) 
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4jte  S  ( l- 1  '+u-u '  ) 


If  |u-u'  | 

7'b 


>  *— J  |  b  > 

5i..  i2  i=i 


h  A*(i'b'  )A+(u'  )A(ib)A(u) 


(III-4) 


where  A^(_/b)  and  k(fo )  have  the  same  meanings  as  A^(K'a')  and  A(Ka). 

In  section  III-l,  V.  is  the  interaction  of  the  atoms  with  radiation 

A 


V  =  -  £  -2s^-  £  * A(r  )  =  -  I  —  £  -A(r  )  +  E  —  £  *A(r  ) 

A  .  m  r»  m  r»  —  -ni 


aJm  c  -  -aj 


o  m^c  ^cr - a  ^-^mc  - oj 


where  A(r)  is  given  by  (II-7).  m  and  m,  p  and  p  ,  and  r  and  r  =r  +p  , 
-  o  a  J‘oJ  ~o  “o  j  “a  ~pj 

are,  respectively,  the  masses,  the  momenta,  and  the  positions  of  the  nucleus 
and  the  J-th  atomic  electron  in  the  o-th  atom.  Neglecting  the  first  sum 
as  compared  to  the  second  one  because  m  »m,  V  can  be  written,  in  second 

C  A. 

quantization,  as 


V  =  Z  <K'a'|—  T  £  •A(r.)|Ka>A+(K’a’)A(Ka), 

A  “  —  —  i  —  —  — 


me  - j 


Assuming  that  the  center-of-mass  coordinates  of  an  atom  coincide  its  nucleus 
coordinates,  the  substitution  of  eqs.  (II-7)  and  (B-5)  gives 


v  =  z  — 

A  me 

Ka 


"jp"  '  BK(K-K’-k){c^(k)ex(k)+C^(-k)e^(-k;}  *<a  '  |  E  e  '-^£  |a> 

1/00  j  =  l  J 


H  A  (K'a' )A(Ka) 


"  1.  k  *  p 

Using  dipole  moment  approximation  e  — j~l,  one  obtains,  after  substituting 


i  r„int  l 
'i  m  tHA  > 

int 

where  H  is  the  hamiltonian  of  the  internal  motion  of  an  atom, 
A  * 


VA  =  ^  "T  *[a+(k)€  (k)+rt(-k)e  (-k)] 

A  yk  c  U  K - a  a  a  a  \ - k  —  \ - X.  — 

Ka  1  LC° 


h  A  ( K ' a ' )A(Ka) 


(HI-5) 


where  ,  =  E  ,-E  is  the  energy  difference  between  the  internal  states 

a  a  a  a 

z 

) a  '  >  and  !a>  and  ed  ,  =  <a  '  |  Z  ep.|a>  is  the  dipole  moment  transition  of 

a  a  j=l  J 

the  atom  from  the  state  !a>  to  the  state  I  a '  >  , 

In  the  same  way,  the  interaction  of  the  free  electron  with  the  radiation 

field 


V  =  —  £.'A(r  ) 

e  J  =  1  me  - J 


can  be  written,  in  second  quantization,  as 


V  =  Z  —  -"C-  (a+(k)e  (k)+a  (-k)e  (-k)).u  5  (u-u'-k)  Af(u')A(u) 

e  \k  mc\  A,  X--K-  “  K - - 


(m-6) 


The  interaction  between  the  ions  and  the  radiation  is  given,  in  second 
quantization,  as 


K£i  \  La 
ib 

7-b-  (111-7 ) 


c 


APPENDIX  C 


NO  PHOTONS  EMITTED  OR  ABSORBED  THROUGH  THE 
INTERACTION  OF  FREELY-MOVING  ELECTRONS  WITH 
A  RADIATION  FIELD 


In  appendix  B,  we  have  obtained  that  the  interaction  of  freely-moving 
electrons  with  a  radiation  field  is  given  by 


eft 


V  =  I 
e  \k  mc  \| 
uu'  * 


2lT^p-  (0^(k)€x(k)+c^(-k)€^(-k)}'u  6K(u-u’-k)A^(u' )A(u), 


L  cn 


Assume  that  the  photons  can  be  emitted  through  the  interaction,  the  number 

7 

of  photons  with  momentum  ftk  and  polarization  k  emitted  per  unit  time  from  1/ 
will  be,  from  section  III-2, 


\(* >  +  Ve  ^r*>  I'  6  WV' 

Let  the  initial  and  final  states  of  the  electrons  and  photons  be 


|(3t}>  =  |...n(u),  n(u' )...>  |  ...T^(k)...> 


J  P  i  rj+ 1  >  =  |...n(u)-l,  n(u' )  +  1. . .  >  |  .  r^(k)  +  !...>, 


then 


n.  (k)  =  Z  L2ne&.  ,  (u-e  )2(f  (k)+l)  N(u)  &„(u-u'-k)  S(E  -E +*») 
A  “  uu'  2  '  —KK  R.  u  u 

m  o)L’ 


2,2  2 


Z  (u-c  )2(f  (k)+l)  N(u)  5(E  -E  +fios)  . 

d  .  o  ~K  K  —  u-k  u 


u  m  0)L 
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Since  the  energy  conserved  delta  function  is  contained  in  the  above  expression, 
a  (k)  is  not  zero  if 


E  ,  -E  +4xjd  =  o 
u-k  u 


i.e. 


y 


cose  =  -  +  7  -  (c-i) 

v  2  u 

where  0  is  the  angle  between  k  and  u  and  v  is  the  speed  of  the  incident  electron. 

q 

Since  — >  1  and  k  and  u  are  positive,  (C-l)  cannot  hold.  Therefore  no  photons 

will  be  emitted  through  the  interaction  of  freely-moving  electrons  with  a  ra¬ 
diation  field.  In  the  similar  way,  one  can  obtain  that  no  photons  will  be 
absorbed  through  the  interaction  of  freely-moving  electrons  with  a  radiation 


field. 


APPENDIX  D 


DERIVATION  OF  EQS .  (lV-31a)  THROUGH  (lV-3If) 

For  derivation  of  eqs .  (lV-31a)  through  (lV-31f),  the  values  of  the 

( *i2) 

following  3j  symbols  taken  from  the  literature w  'are  needed. 


h 

j2 

b 

ml 

nip 

m3 

\m1  m2  my 

1 

0 

i 

0 

0 

0 

*  1/3 

2 

0 

2 

0 

0 

0 

1/5 

2 

1 

1 

0 

0 

0 

2/15 

2 

2 

2 

0 

0 

0 

*  2/35 

3 

1 

2 

0 

0 

0 

*  3/35 

k 

2 

2 

0 

0 

0 

2/35 

0 

1 

1 

0 

1 

-1 

1/3 

2 

1 

1 

0 

1 

-1 

1/30 

0 

2 

2 

0 

1 

-1 

*  1/5 

0 

2 

2 

0 

2 

-2 

1/5 

2 

2 

2 

0 

1 

-1 

1/70 

2 

2 

2 

0 

2 

-2 

2/35 

U 

2 

2 

0 

1 

-1 

8/315 

u 

2 

2 

0 

2 

-2 

_ 1/630 

when  one  needs  the  values  of 


other  than  its  square,  the  negative 


square  root  of  the  number  should  be  taken  if  it  is  preceded  with  the  star 
symbol  *.  Since  the  sum  of  jg  and  for  each  row  in  the  above  table  is 
even,  each  of  the  3j  symbols  in  th  above  is  invariant  in  a  permutation  of 
any  two  columns. 

For  convenience,  we  rewrite  here  eqs.  (IV-20),  (lV-22),  (lV-23),  and 


(IV- 26) 


» 
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119 


G(n,n,q)  =  n^-2Rp  E  <nlm|e^  “|nlm>  +  E  |<nl'm'  |e  ^  “ |nlm>  | 


WE.  i 


I  m 


Im 
I'm' 


(IV- 20) 


G(n,n  Jq)  =  Z  |<n'/'m'|e  ^  —  |  nlm>  | 
Im 
I'm' 


(IV- 26) 

2 


<n'l'm' |e1£  ~|nlm>  =  ,Z  ,Y  ,  ,,  ,R  ,  lfI(q) 

p=  |l-i'|  P,lm,l'm'  p,nl,n  I 


t+r 


(IV- 22) 


Y  ,  ,  p+2tn'  (2p+i)  j  (2l+l)(2l  ,+l)  (VU')  (Vll>)  (IV- 25) 

p,lm,l  m  =  i  v  1  \om-ny  \ooo/ 

where  p  takes  the  integers  between  ll-l'l  and  l+l'  such  that  p+l+l'  =  even, 

otherwise  (  *  Jvanishes. 

\ooo  J 

For  n'=n=l,  it  is  easily  to  obtain  that 


G(l,l,q)  =  {1-R0#10> 


10 


(q)}‘ 


For  n'/n,  one  obtains  from  eq.  (IV- 26)  that 


G(l,2,q)  -  E  Rp^  215+1  ^ 21 ' +1)  (Lo) 


(IV-31«) 


=R0,10,20^  +  5R1,10,21^ 


G(l,5,q)  -  Z  Rt2  .n  ,-|,(q)(2p+l)(2l'+l) 

I'p  P>1 


(v°i\ 

\ooo/ 


(lV-31b) 
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■  E0,10,30('!)  +  3ilO,31(<l)  +  5E2,10,32(,5)  (IV-51C) 


0(2,5, q)  -  ?  R  OA  ,.,(q)(2p+l)(2<'+l) 


i'p  P>20,3i 


/porV 
V  coo  y 


+  E  5R 

r P  p»21 

=  R0, 20,30^)  +  5Ri>20,3l(<l)  +  5r|, 20,32^ 

+  5Rl,21.,30^q)  +  5R0,21,31^q)  +  6R2,21,31^q) 


+  <2l,32(q)+ <21,32^ 


(!V-31e) 


For  n'=n=2,3>  one  obtains  from  eqs,  (IV-20)  and  (IV-26)  that 


G(2,2,q)  =  4-2Re{<200|el£ -|200>  +  <210| e1^  -| 210>  +  2  <21l| e1^ -| 211>  ) 

+  E,  |Rr,  pi  o # '  I  ^ 2p+l)(2l+l)(2f  ,+l)f 
If 'p  P>2I,2I  \ooo  / 

G(3,3,q)  -  9-2Re(<300|  e^  —|300>  +  <310|  ei^'-|310>  +  <320|  e^ -|320> 

+  2  <31l|ei£ -|311>  +  2  <32l|ei£-|32l>  +  2  <322 1 e1-3*’- 1 322 >  ) 

t4lEP,3i,3i'(,5>|2(2ptl)(2m)(2',tl)(coo')2 
The  terms  in  the  brackets  of  the  above  two  equations  are  computed  through 
the  use  of  eqs.  (IV-22)  and  (IV-23).  Then 


G(2,2,q) 


4"2^R0,20,20^ 


2R0, 21,21^ 


2R2,21,21^ 


+  R0,21,21^"2R2, 21,21^^ 


+  E  R«  PO  p..(q)(2p+1){2l'+1)  r 
I  ip  p, 20,22  'ooo 


(p0,')! 


+  E5Rr  Ol  5,.(q)(2p+D(2i 

i'p 


■u)  ry 

\ooo  / 


!'”2R0,20,20^  +  R0,20,20^  "  ^R0,21,21^  +  ^R0,2l,21^ 


+  6Rl,20,21(q)  +  6R2,21,21(q) 


(IV-51J ) 


0(5, 5,q)  -  9-2(R0j}0(30(q)  ♦  3H0j31;31(q>  -  3^>J8>„(q)) 

2 

*(Fp<3^''<,)(2ptl)(2,'n>(-')2 


+  ^  5R  p  , .,(q)(2p+l)(2/’+l) 

|'p 


fp2,,y 

\ooo  y 


9'2R0,30,30(q)  +  R0,30,50(q)  "  6Ro,3l,5i^q)  +  3R0,3l,5l(q) 


10R0,32,32(q)  +  5R0,32,32(q)  +  6Rl,30,3l^  +  12Rl,31,32^q) 


+  6Ro  xi  xt^)  +  ^  Ro  xo  xo(0  +  10R2,3o,32(q)  +  l8R*  21  —(q) 


2,31,31  7  2,32,32' 


3,31,32' 


+  ?R?,32,32(q)- 


(IV-3lf ) 
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IS-  ABSTRACT 

-  3  The  Index  of  refraction  and  the  attenuation  coefficient  per  unit  time,  as 
veil  as  per  unit  length,  has  been  obtained  for  an  arbitrary  medium  using 
Maxwell's  wave  theory.  Results  are  expressed  In  tirms  of  the  microscopic 
currents  due  to  the  motion  of  the  particles  In  the  medium.  The  contribution 
of  the  neutral  atone  to  the  absorption  of  photons  In  plasmas  has  also  been 
Investigated  by  extending  the  Akcasu  and  Wald's  work  on  the  absorption  due  to 
the  Inverse  bremsstrahlung  of  electrons  In  the  field  of  neutral  atoms  to 
higher  electron  temperatures  -and  higher  photon  energies.  This  problem  vas 
formulated  by  using  plane  vaves  for  electron  vave  function.  In  this  vay  an 
expression  has  been  obtained  for  the  absorption  of  coefficient  per  unit 
length  due  to  the  above  absorption  mechanism  In  terms  of  the  elastic  and 
Inelastic  electron-atom  scattering  cross  sections  alloving  the  atoms 
Initially  to  be  In  any  excited  state.  The  absorption  coefficient  has  been  cal¬ 
culated  explicitly  for  hydrogen  atoms;  the  results  are  presented  graphically  as 
a  function  of  electron  temperature  and  radiation  frequency.  (  )  -*£ — 
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